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Abstract 

This set of lectures is an introduction to black-hole astrophysics. The emphasis is 
£Sj , made on the phenomenology of X-ray binaries and of supermassive compact objects at 

galactic centers. 
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1 Introduction 

Black holes (BH) are the most fascinating objects predicted by general relativity. There 
exist about 20 confirmed candidates [I] for astrophysical BHs in the mass range 5-20 M Q 
and about three dozen supermassive BH candidates [2] in the mass range 10 6 - 10 9 ' 5 M & . 

Unfortunately, there exists as yet no direct evidence for astrophysical BHs 0. At present 
we only hope that the black-hole paradigm may be proved or ruled out by comparing BH 
candidates with credible alternatives. Fortunately, BHs are dark and compact, which narrows 
the list of possible alternatives among standard astrophysical objects. For a stellar-mass 
BH, the only standard astrophysical alternative is neutron stars (NS). It is therefore quite 
important to understand the properties of NSs and their observational distinction to BHs. 
Therefore, considerable attention will be devoted here to NSs. 

These lecture notes are in large part based on the standard text books Misner, Thorne, 
and Wheeler [3], Shapiro and Teukolsky |4J, Wald [5], and Carroll and Ostlie [6], and on 
the review articles by Townsend [7j and Narayan [8j . The exception is section [9] which is 
based on original articles. Efforts are made to provide citations to original papers wherever 
appropriate but the list of references is by no means complete. A number of interesting talks 
may be found at the site of a recent conference on supermassive black holes [H] . 

We use the positive-time negative-space signature convention, i.e., (+, ) and we mostly 

use the so-called natural units in which c = % — G — 1. In these units, the physical quantities 
are expressed in powers of the Planck mass mpi = ^/hc/G, Planck length Z P i = ^HG/c 3 , or 
Planck time tp\ = lp\/c. 



2 Preliminaries 

2.1 Spherical Configurations 

Consider the space time metric 

ds 2 = g lxv dx tM dx v , (2.1) 

in which the metric tensor g^ is time independent. In general relativity we distingush 
between static and stationary metric. Stationary configurations are described by the metric 
coefficients that do not depend on time. For static configurations we also require the time 
reversal (t — > —t) invariance of the metric. In this case, the mixed components g^ must 
vanish. 

1 See A. Miiller's lecture at this School 
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Static fluid configurations are spherical. The most general metric generated by a spherical 
mass distribution is of the form 

ds 2 = £ 2 dt 2 - X 2 dr 2 - r 2 (d6 2 + sin 2 Odtf), (2.2) 

where £ and A are functions of r only. The function £, called the "lapse function", may be 
represented in terms of the gravitational potential 

£ = e*« (2.3) 

and the function A is related to the enclosed mass A4(r) 

It may be easily shown [TO] that Einstein's equations reduce to two nontrivial independent 
equations 

2 ^ -_ 87rrA 2 T r r + ^A (2 . 5) 



£ <ir r 

= 87rrA 2 T° - ^ 
A dr r 



l^ = 8nr\ 2 T%-^±. (2.6) 



For a perfect fluid 

V r = -p, T° =p. (2.7) 
Using (|2.4|) . Einstein's field equations take the form 

rf£ = M + 4nr 3 p 

dr * r(r - ' 1 ' J 

4vrr 2 p- (2.9) 



dr 

The latter may be written in the form 



rr 

M(r) = / dr' Anr' 2 p(r') , (2.10) 
Jo 



which shows that the function M. may indeed be interpreted as an enclosed mass. 
Finally, if we impose the particle number conservation constraint ( 1C.1I) . we have 

R 

dr 47rr 2 (l - 2M/r)~ 1/2 n(r) = N, (2.11) 
where we have employed the spherical symmetry to replace the proper volume integral as 

u^dJ:^ = / dr4vrr 2 A. (2.12) 
Jo 
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2.2 Schwarzschild Solution 

Assume the absence of matter in the region r > R. In this region equations ( 12.8ft and ( 12.91) 
may be easily solved. One finds 



M(r) = M = const., (2.13) 



The metric takes the form 

ds 2 = ^ 1 _^ tft2 _ dr 2 -r 2 (d0 2 + sin 2 0# 2 ). (2.15) 

This is known as the Schwarzschild metric. The constant M is the mass of the source. In 
the weak field limit r 3> M, we obtain the Newtonian potential 

M 

</? = ln£^ . (2.16) 

r 

This metric describes the gravitational field outside of any spherical object of mass M, 
including a black hole. The sphere with radius r = 2M, at which the metric diverges, is the 
BH event horizon. 

2.2.1 Birkhoff's Theorem 

Theorem 1 The exterior spacetime of all spherical gravitating bodies (not necessarily static) 
is described by the Schwarzschild metric. 

The proof is simple. See, e.g., Misner et al. [3]. 
2.3 Spherical Stars 

Generally, we assume that p and p satisfy an equation of state, e.g., in the form given by 
(IC.25P and (10.261) . and are generally functions of local T and p which in turn depend on £ 



through the Tolman equations (IC.lOj) . Numerical integration of (12.81) and (12.91) from r = 
to some radius R with the boundary conditions 

f(fl) .( 1 _^) 1/2 i M{ti)=0 ,2,7) 

is rather straightforward and, as a result, a nontrivial spherical distribution of matter is 
obtained. The boundary is usually naturally provided at the radius R where p = p = 0. In 
the case when the p = p = point does not exist, one must integrate up to infinity or up to 
a chosen cutoff radius. 

It is often more convenient (e.g., if the equation of state is given in the form p = p(p)) 
to express the field equations in the Tolman- Oppenheimer-Volkoff (TOV) form 

dp , . M. + 47rr 3 n 

-T = -fr + P)—, 2.18 

dr r(r - 2M) 
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dM n 

— — = 4vrr 2 p. 2.19 
dr 

Here, equation (12.181) is obtained from (12.81) with the help of the equation of hydrostatic 
equilibrium (1B.14|) which may be written as 



din £ 1 dp , 
— j— — -f. 2.20 

dr p + pdr 

The set of equations (I2.18I) - (I2.20I) is called the TOV equations. 
2.4 Newtonian Limit 

General relativity reduces to Newtonian theory in the limit of weak gravity and low velocities. 
The Newtonian limit is achieved by the approximation £ = ~ 1 + M/r < 1 , and 
p ^ p. In this limit, the two equations, ( 12.181) and ( 12.191) . can be combined to give one 
2nd-order differential equation 

1 d r 2 dp , 

/ = -47rp. 2.21 

r A dr p dr 

In this approximation, the relativistic energy density p is just the rest mass density 



p = ran. (2.22) 



2.4.1 Polytropes 



If the equation of state is in the polytropic form (1Q.32j) . the equilibrium configurations are 
called polytropes. Owing to (I2.22p we can rewrite the equation of state as 

p = Kp r , (2.23) 

where 

K = . (2.24) 

It is convenient to write 

T = 1 + l/n, (2.25) 

where n is called the polytropic index. Then the Newtonian equation (I2.2ip can be reduced 
to a simple form by writing 

p = p c 9 ; r = az, (2.26) 

1 1/2 

(2.27) 



a: 



{n + l)Kp, 
Air 



l/n-l 



where p c is the central density. Then 



1 d d6 

^i~ z V = -° ■ 2 - 28 

z z dz dz 
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This is the Lane-Emden equation for the structure of a polytropic index n. This equation 
can be numerically integrated starting from the center with the initial conditions 



9(0) = 1; 0'(O)=O. 



(2.29) 



For n < 5 (T > 6/5), the solutions decrease monotonically and have a zero at a finite value 
z = Z\\ 8(zi) = 0. This point corresponds to the surface of the star, where p = p = 0. The 
radius is 

(n+ l)i^ 1/2 



R = az± 



Air 



pt}- n)/2n zi 



(2.30) 



while the mass is 



M 



4:Trr 2 pdr 



Ana 3 p c 



z 2 9 n dz 



d 2 d6 
—4ira p c / —z —dz 

o dz dz 

3/2 

(3-n)/2n Jla/ I 

p c z iril 



47T 



(n+l)iif 



47T 



(2.31) 



where the absolute value \9[ \ has been introduced because 6[ = 0'(^i) is negative. Eliminating 
p c between (12.301) and (I2.3ip gives the mass-radius relation for polytropes [1]: 



M = 4 n R^-n) /(l-n) 



{n + l)K 

47T 



n/(n— 1) 



(3-n)/(l-n) 21/1/ 



The solutions we are particularly interested in are 

T = 5/3, n = 3/2, Zl = 3.65375, z 2 \9[\ = 2.71406 , 

T = 4/3, n = 3, z x = 6.89685, ^|^| = 2.01824 , 

where the values 5/3 and 4/3 of Y correspond to the nonrelativistic and ultrarelativistic 
regimes, respectively (see appendices IC.4l and 10.51) . 



(2.32) 

(2.33) 
(2.34) 



3 Compact Astrophysical Objects 

Traditionally, compact astrophysical objects represent the final stages of stellar evolution: 
white dwarfs (WD), neutron stars (NS), and black holes (BH). They differ from normal stars 
in two basic ways. 

First, since they do not burn nuclear fuel, they cannot support themselves against grav- 
itational collapse by generating thermal pressure. Instead, either they are prevented from 
collapsing by the degeneracy pressure (WDs and NSs) or they are completely collapsed 
(BHs). With the exception of the spontaneously radiating "mini" BHs with masses less than 
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10 12 kg and radii smaller than a fermi, all these objects are essentially static over the lifetime 
of the universe. 

The second characteristic distinguishing compact objects from normal stars is their ex- 
ceedingly small size. Relative to normal stars of comparable mass, compact objects have 
much smaller radii and hence, much stronger surface gravity. 





Mass 


Radius 


Mean Density 


Surface Potential 


Object 


M 


R 


g cm -3 


M/R 


Sun 


M & 


Rq 


1 


KT 6 


White dwarf 


< M 




< 10 7 


~ 10" 4 


Neutron star 


~ 1 - 3M Q 


~ io- 5 i? 


< 10 15 


~ HT 1 


Black hole 


Arbitrary 


2M 


~ M/R 3 


~ 1 



M Q = 1.989 x 10 30 kg; R Q = 6.9599 x 10 5 km 



3.1 White Dwarfs 

White dwarfs are stars that no longer burn their nuclear fuel and their gravitational collapse 
is supported by the pressure of degenerate electrons. Their mass and radius are about 1 M & 
and 5000 km, respectively. 



3.1.1 Equation of State 

We assume that the interior of a WD is almost completely ionized plasma at a temperature 
T <C m e . Hence, the electrons are assumed to be degenerate. The density p of a WD is 
basically the density of barionic matter (neutrons end protons). The pressure is dominated 
by the pressure of electrons. This may be seen as follows. In the nonrelativistic regime 
the pressure roughly equals the average kinetic energy (see equation (IC.14I) ). Owing to the 
momentum conservation, the electron and proton average momenta are equal and hence, 
their average kinetic energies in the nonrelativistic regime are inverse proportional to their 
masses. Hence, the pressure of electrons is roughly by a factor of m n /m e ~ 2000 larger than 
the pressure due to nucleons. 

To determine the equation of state, first assume that the electrons are nonrelativistic. 
The density is related to the number density of electrons as 

p = rj e m n n e , (3.1) 

where the factor rj e , to a good approximation, equals the number of nucleons per free electron, 
i.e., r] e = A/Z. For example, for completely ionized pure 12 C, r] e = 2. The number of electrons 
is given by ()C.24|) . with m = m e . Hence, 

1 

P = —zr) e m B {m e X) 3 , (3.2) 

where m e X = qp = v/// 2 — ml is the Fermi momentum of the electrons. We could now use 
the relativistic expression (IC.26I) for the pressure together with (13. 2p and numerically solve 
the TOV equations to find the entire range of white-dwarf solutions. 
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However, it is instructive to consider the nonrelativistic and extreme relativistic regimes 
separately because, in these cases, the equation of state takes a polytropic form. The pressure 
in the two regimes is given by ( 12.231) with T = 5/3 for the nonrelativistic and T = 4/3 for 
the extreme relativistic regime. Using this and the solutions for the polytropes f l2.30l) -( l2.34l) 
we find the radius and the mass of the WD: 



Low-density (nonrelativistic) regime. 



r 4 



R = 1.122 x 10 4 



K 



3 2 / 3 7T 4 / 3 
c 5 / 3 5 / 3 

5m e m n rje 



Pc 



10 6 g cm - 



-1/6 



km, 



(3.3) 
(3.4) 



M = 0.4964 



0.7011 



10 6 g cm 

(—)' 

Vl0 4 kmy 



1/2 



-5/2 



M, 



© ) 



M 



© 



(3.5) 



High-density (ultrarelativistic) regime. 



R = 3.347 x 10 4 



K 



Pc 



3 1 / 3 7T 2 / 3 



10 6 g cm 



-1/3 



*r 2/3 km, 



M 



1.457 [-^J M . 



(3.6) 

(3.7) 
(3.8) 



Note that M is independent of p c in the extreme relativistic limit. We conclude that 
as p c — » oo, the electrons become more and more relativistic and the mass asymptotically 
approaches the value (I3.8P as R — » 0. This mass limit is called the Chandrasekhar limit, and 
represents the maximum possible mass of a WD [11] . 

Of course, the limiting radius is not zero. The integration of the TOV equations using 
the exact degenerate Fermi gas equation of state would also give the limiting value of R of 
the order of 

#ch ~ ~ 5 x 10 3 km. (3.9) 

m n m P 



3.2 Neutron Stars 

If the mass of the collapsing star is larger than the Chandrasekhar limit, the degeneracy 
pressure of the electrons can no longer support the gravitational attraction and the collapse 
does not stop. As the density increases, the Fermi energy E-p = p of the electrons increases 
according to (13. 2p . 



9 



At a density of about 2 x 10 7 g cm -3 , the Fermi energy of the electrons has risen to 
m n — m p = 1.29 MeV where electrons can now be absorbed by protons through the inverse 
(3 decay 

e~ + p — > n + v e (3.10) 
This reaction cannot come to equilibrium with the reverse reaction 

n + v e — > e~ + p (3-11) 

because the neutrinos escape from the star and the normal (3 decay cannot occur because all 
electron energy levels below E = m n — m p are occupied when Ep > m n — m p . At a density 
in the range 10 7 < p < 4 x 10 11 g cm -3 the medium is a composition of separated nuclei in 
equilibrium with a relativistic electron gas. At p ~ 4 x 10 11 the ratio n/p reaches a critical 
level. Any further increase leads to a "neutron drip" - that is , a two-phase system in which 
electrons, nuclei, and free neutrons coexist. When the density exceeds about 4 x 10 12 g cm -3 , 
more pressure is provided by neutrons than by electrons. The neutron gas so controls the 
situation; one can describe the medium as one vast nucleus with lower-than-normal nuclear 
density. 

As the density reaches the normal nuclear density of about 

Pnuci = t^-^ = 2 x 10 14 gcm- 3 , (3.12) 

there is a phase transition in which nuclei dissolve. The resulting fluid consists mostly of 
neutrons with a small ( ^ 5%) fraction of electrons and protons. The pressure is dominated 
by the degenerate (T = 0,/i n > m n ) gas of neutrons. Can the degeneracy pressure of 
neutrons support the star against collapse? 

Assuming that the Newtonian calculations are still valid at a density not very much 
higher than p nuc i, we can use the high-density results obtained for WDs but with m e m n . 
The maximal mass remains the same as in (13.91) but the value of the critical R is of the order 
of 

R „ ^ ~ 2.5km, (3.13) 

mi 

which is close to the Schwarzschild radius of the sun, so the neglect of GR effects is not 
justified. 

We have to solve the TOV equations with the equation of states of nuclear matter. The 
first numerical calculation was performed by Oppenheimer and Volkoff (12] who used the 
relativistic degenerate Fermi gas described by flC.24p - flC.26p . In Fig. [T]we plot the mass of 
the NS as a function of the radius R. The maximum of the curve corresponds to the OV 
limit. The limiting values are [13] 



^ = 3.357 (!) " 2 = 9.6 2 (hm) ^ , (3,4) 

m \9j \ m 2 J \g J 



in 



3 '2\ 1/2 n ( 1 GeV\ 2 /2\ 1/2 



M ov = 0.38426 I- =0.7 I— ( - 1 M , (3.15) 



10 



M 



0.40- 



o . i o -; 



30 - 



20 - 




2 



4 



6 



8 



R 



Figure 1: Mass versus radius for fermion stars at zero temperature in the general- relativistic 
framework (solid line) compared with the corresponding Newtonian approximation (dotted 
line). The dashed line is the Schwarzschild BH limit M = R/2. 

where g is the fermion degeneracy factor. This limit is reached when the central density 
becomes p c = 5 x 10 15 g cm -3 with g = 2 for neutrons. The part of the curve left from the 
maximum represents unstable configurations that curl up around the point corresponding to 
the infinite central-density limit. 

However, the degenerate neutron gas equation of state is not realistic because at such 
large densities the effects of nuclear forces must be included. More realistic equations of 
state predict the maximum NS mass in the range 1.5 - 2.7 M & . Hence, the maximal mass 
is rather sensitive to the not very well-known equation of state for nuclear matter. Since 
the density inside the star varies from very large central values to zero at the surface, the 
equation of state is actually rather complicated as the star may contain different phases of 
nuclear matter. 

The possibility to identify some compact objects as BHs relies in part on being able to 
state categorically that the observed object has a mass larger than the maximum allowed 
mass of a stable NS. It turns out that it is possible to set an upper limit to the NS mass 
based on rather general assumptions (Rhodes and Ruffini [Tl]): 

1. The TOV equation determines the equilibrium structure. 

2. The equation of state satisfies the local stability condition 



that is, the speed of sound is real. If this condition were violated, small elements of 
matter would spontaneously collapse. 




(3.16) 
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3. The equation of state satisfies the causality condition 

c 2 s < 1, (3.17) 
that is, the speed of sound is less than the speed of light. 

4. The equation of state below some "matching density" po is known. 

Rhodes and Ruffini performed a variational calculation to determine which equation of 
state above po maximizes the mass. Then, the numerical integration of the TOV equations 
for a chosen equation of state below p gives 

/ \ -i/2 

M max ~3.2 P l M . (3.18) 

\4.6 x lO^gcm - ' 3 / 

A semianalytic treatment of Nauenberg and Chapline fl5j with similar assumptions about 
the equation of state gives 

M max ~3.6M . (3.19) 

Abandoning the causality constraint still leads to a severe mass limit, assuming general 
relativity to be valid. One finds (see, e.g., Hartle and Sabbadini |16j) 

M max ~ 5.2M . (3.20) 

In their estimate, Rhodes and Ruffini took the the so-called Harrison- Wheeler equation 
of state which accurately describes the nuclear matter densities below the neutron drip. In 
fact, it turns out that the upper mass limit is not very sensitive to the equation of state used 
below p . It is important to note that even if a new physics exists at subnuclear and subquark 
level (preons, pre-preons etc.; see, e.g., a black-hole sceptical paper [H]), it is reasonable to 
expect that the equation of state will still satisfy the above conditions and hence, the above 
limits cannot be significantly altered. 

3.3 Black Holes 

What happens if the mass of the collapsing star is larger than the maximal allowed NS mass? 
In such a case, there is nothing to prevent the star from further collapse ending in a BH0. 

Astrophysical BHs are macroscopic objects with masses ranging from several M (X-ray 
binaries) to 10 6 — 1O 9 ' 5 M (in galactic nuclei). Being so massive, these BHs are described 
completely by classical general relativity. As such, each BH is characterized by just three 
numbers: mass M, spin parameter a, defined such that the angular momentum of the BH 
is J = aM, and electric charge Q [3J HI EH] (see also J. Zanelli's lectures at this School). 
Actually, an astrophysical BH is not likely to have any significant electric charge because it 
will usually be rapidly neutralized by surrounding plasma. Therefore, the BH can be fully 
characterized by measuring just two parameters, M and a, of which the latter is constrained 
to lie in the range from (nonrotating BH) to M (maximally-rotating BH). 

2 It has been recently proposed that loop quantum gravity effects stop the collapse to a singularity by a 
bounce of the infalling matter [TS] 
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3.3.1 Spherical Collapse 

A useful toy model that illustrates the collapse is a self-gravitating spherically symmetric 
ball of dust (i.e., zero pressure fluid). Birkhoff's theorem (section 12.2. ip implies that the 
metric outside the star is Schwarzschild (section 12.21 equation ( 12. 15ft ). This is valid outside 
the star but also, by continuity of the metric, on the surface. If r = R(t) on the surface, we 
have 



ds 2 



2M 
~R 



2M 
~R 



-i 



fdR\' 



dt 2 - R 2 (d6 2 + sin 2 9dcf> 2 ). 



(3.21) 



Zero pressure and spherical symmetry imply that a point on the surface follows a radial 
timelike geodesic, d9 = d(j) = and ds 2 = dr 2 > 0, so 



1 - 



2M 



2M 



-i 



fdR 
\~dt 



i 2 



where 



i- dt 
dr 



(3.22) 



(3.23) 



Since £ = d/dt is a Killing vector, by Proposition [T] in appendix IA.5I the energy per unit 
mass 



E = i, 



dx" 
dr 



9oot 



1 



2M 



t. 



(3.24) 



is a constant of motion. Note that E > for timelike and null geodesies as long as £ M is 
timelike, i.e., as long as R > 2M. The quantity E is constant along the geodesies and E < 1 
for gravitationally bound particles. Using (I3.24[) in (I3.22|) gives 



fdR 
\~dl 



2M 
~R 



2M 
~R 



1 + E' 



(3.25) 



We plot this function in Fig. [2j The surface radial velocity dRj dt as a function of R has a zero 
at R = R max and a minimum at R = 2M. We consider the collapse to begin at R = -R max 
with zero velocity. The radius R then decreases and approaches R = 2M asymptotically as 
t — > oo. This may be seen by integrating 



t = E 



dR(l - 2M/R)-\2M/R - 1 + E 2 )~ 1/2 , (3.26) 
2M. So an observer "sees" the star contract 



where -R m i n > 2M. Obviously t — > oo as -R m i n 
at most to R = 2M but no further. 

However, from the point of view of an observer on the surface of the star, the relevant 
time variable is the proper time along a radial geodesic, so we use 



d _ 1 d If 2M 
~db ~ tdr ~ E V ~ 



d 

dr 



(3.27) 
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Figure 2: The surface radial velocity squared versus R for a spherically symmetric collapse 
of a ball of dust (figure from [7]). 




Figure 3: The surface radial velocity as seen by the observer on the surface (figure from 
0)- 



14 



to rewrite f!3.25j) as 



The star surface falls from R = -R max through R = 2M in finite proper time. In fact, it 
reaches the R = singularity in the proper time 

So, for an observer following radial geodesies, the star collapses to a singularity in a finite 
proper time of the order of 10 -5 s (~ Rsch/c — 3km/300 000 km s _1 ). 

Nothing special happens at R = 2M, which suggests that we investigate the spacetime 
near R = 2M in coordinates adapted to infalling observers. 



3.3.2 Eddington-Finkelstein Coordinates 

A singularity of the metric is a point at which the determinant of either the metric tensor or of 
its inverse vanishes, or at which some elements of the metric or its inverse diverge. However, 
a singularity of the metric may be simply due to a failure of the coordinate system. A simple 
two-dimensional example is the origin in plane polar coordinates. Such singularities, usually 
referred to as coordinate singularities, are removable by a coordinate transformation. 

If no coordinate system exists for which the singularity is removable, then it is irre- 
movable, i.e., a genuine singularity of spacetime. Any singularity for which some scalar 
constructed from the curvature tensor blows up as it is approached, is irremovable. Such 
singularities are called curvature singularities or true singularities. The singularity at r = 
in the Schwarzschild metric is an example. 

We now show that the apparent singularity of the Schwarzschild metric at the Schwarz- 
schild radius is removable. On radial null geodesies in Schwarzschild spacetime 

it2 = lT^f ir ^ (ir " f ' (3M) 



where 

r-2M 

r* = r + 2M In 



2M 



(3.31) 



is the Regge-Wheeler or tortoise radial coordinate. As r ranges from 2M to oo, r* ranges 
from — oo to oo. Thus 

d(t±r*) = (3.32) 

on radial null geodesies. Now, define a new time coordinate (ingoing radial null coordinate) 
v by 

v = t + r*, — oo < v < oo (3.33) 
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and rewrite the Schwarzschild metric in ingoing Eddington-Finkelstein coordinates |20j v , r, 
9, and (f>: 



ds 2 



l-==-J (dt 2 -dr* 2 ) -rW 
1-— ^ dv 2 -2drdv-r 2 dQ 2 . 




(3.34) 



This metric is initially defined for r > 2M because the relation v = t + r*(r) between u 
and r is defined only for r > 2M, but it can now be analytically continued to all r > 0. 
Because of the dr dv cross-term, the metric in EF coordinates is nonsingular at r = 2M, 
so the singularity in Schwarzschild coordinates is really a coordinate singularity. There is 
nothing at r = 2M to prevent the star from collapsing through r = 2M. 

However, no future-directed timelike or null worldline can reach r > 2M from r < 2M. 
This may be seen as follows. When r < 2M, we have 



Hence, for all timelike or null worldlines, dr dv < 0. The null coordinate dv > for future- 
directed worldlines, so dr < with equality when r = 2M, dQ = (i.e., ingoing radial null 
geodesies at r = 2M). 

Thus, no signal from the star's surface can escape to infinity once the surface has passed 
through r = 2M. The star has collapsed to a black hole. For the external observer, the 
surface never actually reaches r = 2M, but as r — > 2M, the redshift of light leaving the 
surface increases exponentially fast [19] and the star effectively disappears from view within 
a time ~ M/m 2 ^ ~ 5 x 10 -6 s . The late time appearance is dominated by photons escaping 
from the unstable photon orbit at r = 3M. 

From the point of view of an outside observer, a star collapsing to a black hole never 
appears to collapse, but rather freezes at the horizon. How then can it be said that the star 
collapses to a singularity, if it never appears to collapse even till the end of the Universe? 

The star does in fact collapse inside the horizon, even though an outside observer sees the 
star freezing at the horizon. The freezing can be regarded as a light travel time effect. Space 
can be regarded as falling into the black hole, reaching the speed of light at the horizon, 
and exceeding the speed of light inside the horizon [2T]. The photons that are exactly at the 
horizon and are pointed radially upwards stay there for ever, their outward motion through 
space at the speed of light being canceled by the inward flow of space at the speed of light. 
It follows that it takes an infinite time for light to travel from the horizon to the external 
observer. The star does actually collapse: it just takes an infinite time for the information 
that it has collapsed to get to the outside world! 




when ds 2 > . 



(3.35) 
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4 Rotating Black Holes 



The spacetime around a rotating black hole is described by the Kerr metric expressed in 
Boyer-Lindquist coordinates 0, BJ E] 

2Mr\ , 2 AaMr sin 2 9 , , , 
ds A = (1 — dt 2 + dt dcj) 

2-i J Zj 

(r 2 + a 2 ) 2 - Aa 2 sm 2 6\ . 2 , 2 E, 2 



sin 2 0d0 2 - — c/r 2 - Edr, (4.1) 
with 

A = r 2 -2Mr + a 2 , E = r 2 + a 2 cos 2 . (4.2) 

The parameters a and M are related to the total angular momentum J = aM. The horizon 
occurs at those points where A = 0, i.e., at the roots of the quadratic equation A = 



r± = M± VM 2 - a 2 . (4.3) 

Note that \a\ must be less than M for a black hole to exist. If a exceeded M, one would 
have a gravitational field with a "naked" singularity, i.e., one not within an event horizon. 
A naked singularity is forbidden by the so called cosmic censor conjecture (for details, see 
Wald [5 J and references therein). A black hole with \a\ = M is called a maximally rotating 
black hole. 

Remarks 

• When a = 0, the Kerr solution reduces to the Schwarzschild solution. 

• Taking — > —0 effectively changes the sign of a, so we may choose a > without loss 
of generality. 

4.1 Geodesic Motion 

A straightforward approach to finding all geodesies is to integrate the geodesic equation 
(IA.12I) directly. However, it is often more economic to simplify the problem by symmetry 
considerations. The Kerr spacetime is stationary and axially symmetric so there exist two 
Killing vectors £ and if) which by Proposition [1] (appendix IA.5j) yield a conserved energy E 
and an angular momentum L (per unit mass) along geodesies 

/ 2Mr\ . 2Marsm 2 6- , A AS 

E = u^=ll- — U + 0, (4.4) 



2Marsin 2 fl. ( [r 2 + a 2 ) 2 - Aa 2 sin 2 e\ . 2 • 
L = -?AV = >— j sin 2 00, (4.5) 

where x M = w M = dx^/dr. In addition we have 

g^x u = k, (4.6) 
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where k — 1, 0, and — 1 for timelike, null, and spacelike geodesies, respectively. One may 
use equations (I4.4p and H4.5[) to solve for i and <p in terms of E and L, and substituting the 
results into ( 14.6[) one obtains a differential equation for f. In the case of equatorial geodesies, 
= tt/2, one finds 

}-f + V(r) = 0, (4.7) 



where 



M L 2 l r ^ 2 / a 2 \ M 



Thus, the problem of obtaining the geodesies in the equatorial plane reduces to solving a 
problem of ordinary, nonrelativistic, one-dimensional motion in an effective potential. The 
calculations are relatively simple for circular orbits. Circular orbits occur where r = 0, which 
requires 

dV 

V = 0, — = 0. (4.9) 

For k > 0, equations (14.91) have solutions E{r) and L{r) for all r > r ph (for details, see 
[U[22]), where r > r ph is the radius of the photon circular orbit. 

4.1.1 Photon Circular Orbit 

Photon circular orbits are possible only for particular radii r p h- The requirements (14.91) for 
null geodesies (k = in (14.81) ) yield a cubic equation in y/r 



r 2 - 2M ± 2aVWr = (4.10) 

with the solution [22] 

~2 



ph 



2M 1 + cos 



^ cos~ 1 (=Fa/M) 



(4.11) 



For a = 0, r p h = 3M, while for a = M, r p h = M (corotating) or r p h = 4M (counterrotating 
orbit). For k = 0, equations (14.91) have a solution 



±E x 3r 2 + a 2 . (4.12) 



The photon orbit is the innermost boundary of the circular orbits for particles, i.e., for 
timelike geodesies (k > 0). 

4.1.2 Angular Velocity 

The angular velocity of an orbiting particle is defined as 
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Using the parameterization flB.9j) for the velocity components, it may be easily shown (ex- 
ercise) that fl can be expressed as 

n = _ gb» + A gto > (4M) 
9<j>(j> + A9o<t> 



where 



A = -^ = -^- (4.15) 



For equatorial circular orbits one finds 



E u Q 



M 1 ' 2 



This is the general-relativistic form of Kepler's third law for equatorial circular orbits. In 
this case, the quantity fl is called the Keplerian frequency. 

4.1.3 Innermost Stable Circular Orbit 

It is obvious that not all circular orbits will be stable because, in addition to (14.91) . stability 
requires 

d 2 V 

frT " °- (4 - 1?) 

From (14. 8p we obtain 

l-£ 2 >— . (4.18) 

Substituting the solution E(r) for timelike geodesies, we obtain a quartic equation in yfr for 
the limiting case of equality. The solution r is is the radius of the innermost stable circular 
orbit (IS CO) also referred to as the marginally stable orbit [22] 



M [3 + Z 2 t (3 - Z 1 ) 1/2 (3 + Zt + 2Z 2 ) 1/2 } , (4.19) 



with 



Zi = 1 + I 1 - ] 



2 \ 1/3 



a 2 \ 1/3 / a 2 \ 1/3 ' 



M 2 J [V 1 + ^V + V M 2 ) 

a 2 \ 1/2 

3^ + Zi J • (4-20) 



For a = 0, r is = 6M; for a = M, r is = M (corotating) ri S = 9M (counterrotating). Obviously, 
compared with the photon circular orbit, the innermost stable circular orbit satisfies ri S > r p h- 
A quantity of great interest for the potential efficiency of a BH accretion disk as an energy 
source is the binding energy of the ISCO. Defining the efficiency r] as the maximum binding 
energy per unit rest mass, from (I4.18P (with =) one finds 

n = l-E is =l-(l-— J . (4.21) 
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Plugging in the solution (14.191) . one finds that the efficiency 77 increases from 1 - y/879 
(a=0) to 1 — a/1/3 (a=M) for corotating orbits, while it decreases from 1 — a/8/9 (a=0) to 
1 — a/25/27 (a=M) for counterrotating orbits. The maximum binding energy for a maximally 
rotating BH is 42.3% of the rest-mass! This is the amount of energy that is released by 
matter spiraling in toward the BH through a succession of almost circular equatorial orbits. 
A negligible amount of energy is released during the final plunge from ri S into the BH. 



4.2 Ergosphere 

A curious property of rotating BHs is that there exist particle trajectories (i.e., timelike 
geodesies) with negative energies. The energy defined by (14.41) can be negative only if the 
time translation Killing vector £ is spacelike. The vector £ is normally timelike at 00 and 
in Schwarzschild spacetime it is timelike everywhere outside the horizon. However, in Kerr 
spacetime it need not be timelike everywhere outside the horizon because 

/ 2Mr \ 

= 200 = 1 - -J— ~ 2 ~n > ( 4 - 22 ) 

\ r z + a z cos z 9 J 



so £ is timelike only if 
This implies 



r 2 + a 2 cos 2 9 - 2Mr > 0. (4.23) 
r < M - VM 2 - a 2 cos 2 9, 



> M + VM 2 - a 2 cos 2 9. (4.24) 



The inner boundary of this region is not physically relevant as it is beyond the horizon. The 
outer boundary , i.e., the hypersurface 



r = M + VM 2 - a 2 cos 2 9 (4.25) 

is called the ergosphere or the stationary limit surface. The ergosphere intersects the event 
horizon at 9 = 0, tt, but it lies outside the horizon for other values of 9. Thus, £ can become 
spacelike in a region outside the event horizon. This region is called the ergoregion. Hence, 
a particle trajectory inside the ergosphere may have negative energy! 



4.2.1 Penrose Process 

In 1969 Penrose exploited this property of Kerr BHs to design a mechanism for the extraction 
of energy from a BH [23]. The mechanism proposed by Penrose can be understood as follows. 
Suppose that a particle approaches a Kerr BH along a geodesic. If g M is its 4-momentum, 
we can identify the constant of motion 

E = q% (4.26) 

as its energy (see section lA. 5 1) . The trajectory is chosen so that it penetrates the ergosphere. 
Now suppose that the particle decays into two others, one of which falls into the hole with 
energy E in , while the other escapes to 00 with energy E ont . By conservation of energy 

E ont = E — E in . (4.27) 
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Normally, E- m > 0, so E oxA < E, but E- m = is not necessarily positive in the ergoregion 
because £ may be spacelike there. Thus, if the decay takes place in the ergoregion, we may 
have E out > E, so energy has been extracted from the black hole. 

The energy extraction by the Penrose process is limited by the area theorem of BH 
mechanics, which states that the surface area of the BH horizon never decreases [21]. The 
area of the horizon is 

A = J ^detg dOdcj) = 8vrM (m + VM 2 - a 2 ) , (4.28) 

where det g = geeg^ is the determinant of the metric on the horizon surface. The surface 
metric is obtained from ( 14. ip by setting dt = dr = 0, r = r + = M + V 'M 2 — a 2 . The 
maximum energy extracted by the Penrose process is obtained if the BH spin reduces to 
zero, i.e., if the BH becomes Schwarzschild. A Schwarzschild BH with the same area will 
have a mass Mj rr (irreducible mass) which satisfies 

167rMi 2 rr = 8nM (m + VM 2 - a 2 ) . (4.29) 

This gives M 2 r = M/2 for a maximally rotating BH (a = M). Hence, the maximally 
extracted energy from a BH of mass M is M - M irr = M(l - 1/V2). This represents w 29% 
of the mass energy of the BH. 

The Penrose process for extracting energy from a rotating BH requires particular con- 
ditions [22] which are very difficult to realize in nature. A more promising mechanism for 
extracting energy from a rotating BH is via magnetic fields (Znajek and Blandford [25l [26] .) 
because magnetic fields are capable of connecting regions very close to the BH to regions 
farther out. Recent general-relativistic magneto-hydrodynamics (MHD) simulations of a 
magnetized plasma in the ergosphere of a rotating BH (Koide et al. [27J [28j [29] , Semenov et 
al. [30], de Villiers et al. [31], McKinney and Gammie [32]) show that, near the equatorial 
plane, the field lines are azimuthally twisted and the twist then propagates outward and 
transports the energy along the rotation axis. The process is dubbed the MHD Penrose 
process because of a close analogy between this mechanism [2H] and Penrose's original idea 

5 Stellar Mass BHs versus NSs 

Astrophysical observation of NSs and BHs is difficult because they are dark and very com- 
pact. Besides, a BH the mass of which is of the order of a few solar masses is not easily 
distingushed from a NS because of their similar properties. So far, NSs and stellar mass BH 
candidates have been identified only as the so-called X-ray binaries. 

5.1 X-ray Binaries 

An X-ray binary is an X-ray source with an optical companion, usually a normal star. Most 
of the Galactic X-ray sources are probably compact objects accreting gas from the companion 
star. 

This interpretation of the observational data follows from these facts [1]: 
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1. The variability of X-ray emission on short timescales implies a small emitting region. 

2. Many of the sources are positively confirmed to be in binary systems, with optical 
primaries orbiting optically invisible secondaries. 

3. Mass accretion onto a compact object, especially a neutron star or a black hole, is 
an extremely efficient means of converting released gravitational potential energy into 
X-ray radiation. 

In general, the list of possible Galactic X-ray source candidates includes all three kinds 
of compact objects: WDs, NSs, and BHs. But in special cases, the specific nature of the 
compact object can be identified. 

5.1.1 Binary X-ray Pulsars 

Binary X-ray sources displaying periodic variations are called binary X-ray pulsars. The 
pulse periods are observed in the range 0.5 s P % 1000 s. Those with short periods of 
about 1 s are normally identified with rotating NS. 

The standard model explains X-ray emission as due to the conversion of the kinetic energy 
of the accreting matter (coming from the intense stellar wind of the companion optical star) 
into radiation, because of the interactions with the strong magnetic field of the neutron 
star, of the order of 10 7 -10 9 T. The magnetic field of the compact object drives the accreted 
matter onto the magnetic polar caps, and if the magnetic field axis is not aligned with the 
spin axis, then the compact object acts as a "lighthouse", giving rise to pulsed emission 
when the beam (or the beams, according to the geometry) crosses our line of sight. 

The reasons why WD and BHs cannot be pulsars are: 

• Rotating WD are excluded because they are too large and their density too low to 
explain such short periods. 

• Rotating BH are excluded because they are axially symmetric and have no structure 
on which to attach a periodic emitter. Any mechanism depending on accretion would 
not be periodic to the observed precision. However, accretion disks round black holes 
could produce the so-called quasi periodic oscillations (see section 15.31) 

5.1.2 Mass Measurements 

The first important criteria that may be used to distinguish a BH from a NS are their masses. 
The fact that there exists a maximum mass of a compact relativistic star such as a NS of 
the order ~ 3M (see section. 13 . 2f) allows the following simple criteria for identifying BH 
candidates [8]: 

If a compact astrophysical object has the mass larger than about 3M & , then the object is very 
likely a black hole. 

The most reliable means of determining astronomical masses are via Kepler's Third Law. 
Consider two spherical masses Mi and M2 in a circular orbit about their center of mass 

3 Obtained from conservation of the magnetic flux during the process of collapse from a "normal" star 
(B ~ l(T 3 -l(r 2 T, R ~ 10 6 km) to a neutron star (R ~ 10 km) 
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(CM). The separation of the two masses is a and their distances from the CM are a\ and 
a 2 . Clearly, a = a\ + a 2 and M\d\ = M 2 a 2 by the definition of the CM. Any spectra emitted 
from, e.g., M 2 will be Doppler shifted, depending on the orbital velocity projection v 2 of M 2 
along the line of sight: 

v 2 = — — a 2 smt , (5.1) 

-'orb 

where P Q rb is the orbital period and i the inclination of the orbital plane to the line of sight. 
Thus, if the spectrum of M\ shows periodic variations, then P orD and v 2 can be measured and 
hence one gets a 2 sini. Alternatively, for X-ray pulses one can measure periodic variations 
in the time of arrival of pulses. The amplitude of these variations is simply the light travel 
time across the projected orbit - that is, a 2 sinz/c. 
Now, Kepler's Third Law states 

G(M 1 + M 2 ) _ ( 2vr , 



a 3 \Porb 

Note that this is valid also for elliptical orbits in which case a is the semimajor axis of the 
ellipse. Using this and 

M, + M 2 

a = — — a 2 , (5.3) 

we obtain 

f(M) - (Mi sin i) 3 _ P OTh vj 
/(Ml) = {Mi + ~ ~^G' 
The quantity / is called the "mass function" and depends only on the observable quantities 
P or b and v 2 (or a 2 sini). 

For several X-ray binaries, it has been possible to measure the mass functions fo = /(Mo) 
and fx = f(Mx) for both the optical companion and the X-ray source, respectively. The 
ratio fo / fx gives th mass ratio 

Mo //o Nl/3 



and then from (15 .4p we can write 



M x = f x { -^L. (5.6) 
sm 1 

A unique value for Mx still depends on knowing sinz. In practice, the X-ray eclipse duration 
and/or variation in the optical light curve are used to set geometrical constraints on sini. 

In this way a complete description of the binary system has been obtained for six eclipsing 
X-ray pulsars with optical companions [4J. Their masses range from 1 to 2.3 M , with error 
bars such that all the data fit in the range 1.2 - 1.6 M , which is expected on the basis of 
current theoretical scenarios for NS formation. One of these X-ray pulsars is the famous 
Hercules X-l neutron star discovered in 1972 [33] in the data of the first astronomy satellite, 
Uhuru launched by NASA off the coast of Kenya@ 



"Uhuru" means "freedom" in Swahili; the launch occurred on the anniversary of Kenya's independence. 
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5.2 Black-Hole Binaries 



Black hole X-ray binaries, or short, BH binaries are the binary X-ray sources with observed 
masses larger than 3 M & with nonperiodic time variability. As of today there are a total of 
20 confirmed BH binaries (Remillard and McClintock pQ) Their large masses makes them 
strong candidates for BHs. However, the mass estimates are reliable only for those for which 
the inclination angle i is well known, which is not always the case. Fortunately, according 
to (15. 6p . M > f(M). The mass function f(M), which depends only on the two accurately 
measured quantities t>2 and P rb, is a strict lower bound on M. Most of the 20 X-ray binaries 
have f(M) itself larger than or of the order 3M . Therefore, these systems are excellent BH 
candidates, regardless of uncertainties in their inclinations and companion star masses. 



5.2.1 Cygnus X-l: A Black-Hole Candidate 

Cygnus X-l is a typical BH binary discovered in 1972 (Webster and Murdin [31], Bolton 
[35] ) as a first stellar-mass BH candidate. The X-ray source (X) of Cyg X-l is variable on 
all timescales varying from ms to months and years. Recent observations show that Cyg X-l 
periodically cycles through two accretion or spectral states: hard and low X-ray states [36J. 
The most dramatic variability is the 1-ms bursts, which set a maximum size for X of the 
order R ^ 300 km and establish the object to be highly compact. 

The optical companion star (O) in the Cyg X-l system is a typical supergiant star with 
a well-known spectrum and the mass of at least Mo ~ 8.5M [I]. From the measurements 
of the orbital elements, P D rb = 5.6 days and ao sin i = (5.82 ± 0.08) x 10 6 km, one obtains 
the mass function of X, fx = (0.252 ± O.O1O)M . Setting sinz = 1 in (I5.6p . one obtains a 
minimum value for the mass 

M x > 3.3M . (5.7) 

It is possible to set a convincing lower limit assuming nothing about the mass of O and 
using only the absence of a prominent X-ray eclipse and the estimate of the radius of O [37J. 
The absence of eclipse implies 

R 

cosi > — , (5.8) 
a 

where R is the radius of O and a = ax + ao is the separation of the two objects. Using (15. 3j) . 
we can write this as 

«"'> R -^Z ■ (5-9) 
~~ a sin i M x + M Q v ' 

Hence, 

,2 

[ao smz 

The function sin x cos x 2 has a maximum value of 2/ (3 a/3), and so 



f R 2 

M x sini cos 2 i > ^ ; ^ ,, 2 ■ (5.10) 



M x >— ( I 5 - 11 ) 

2{ao simy 



Estimating R from luminosity and the effective temperature, one finds [1] 

' d 
1 kpc 



R = 6.62 x 10 6 ( — — ) km. (5.12) 
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From (I5.12p with known fx and do sin i one finds 



d ^ 2 



M x > 3.4 (gj— J M 0- ( 5 - 13 ) 

The distance c? to Cyg X-l is determined using two methods: a) from the assumed luminosity 
compared with the apparent luminosity of O and b) from the absorption vs. distance curve 
calibrated from a large sample of stars in the same direction. The estimated distance for O is 
d ~ 2.5 kpc with an absolute minimum of 2 kpc in order to produce the observed absorption. 

We can summarize the situation as follows: the lower limit (15 . 1 3[) of 3.4M is very solid. 
Adopting the more reasonable value of d ~ 2.5 kpc increases this to 5.3 M & . Various other 
less rigorous but more realistic arguments, as well as more recent measurements of R, d, and 
fx give even larger values. The currently accepted mass range of Cyg X-l is 6.8-13.3 M 

HI- 



5.3 Spin Estimates 

Measuring the BH spin amounts to measuring the Kerr-spacetime parameter a. In contrast 
to BH mass estimates where Newtonian gravity applies, the spin of a BH or of any other 
rotating astrophysical object does not have any Newtonian effect on the surrounding objects. 
Only for relativistic orbits does spin have measurable effects. Therefore, to measure a, we 
need test particles orbiting very close to the innermost stable circular orbit (ISCO) (see 
section 14.1.31) . Such test particles are provided by the accretion disk. 

The gas in an accretion disk starts from large radii and spirals in through a sequence of 
nearly circular orbits as it viscously loses angular momentum. The main source of instability 
and loss of angular momentum are MHD effects, e.g., magnetorotational instability [38J and 
magnetoviscous instability |39j. When the gas reaches the ISCO, no more stable circular 
orbits are available, so the gas accelerates radially and free-falls into the BH. Thus, the ISCO 
serves effectively as the inner edg<S of the accretion disk. A variety of observational methods 
have been proposed for estimating the radius r in (spectral fitting method; relativistic iron 
line method), the Keplerian frequency Q (quasiperiodic oscillations method), or the binding 
energy 77 at r- m . Next, we briefly describe some of these methods. For more details and for 
references, see [8]. 



Spectral Fitting 

When a BH has a large mass accretion rate M, corresponding to an accretion luminosity 
L acc above a few per cent of LEdd (see section 16.3.11) . the accreting gas tends to radiate 
approximately as a blackbody. In this spectral state one can theoretically calculate the flux 
of radiation F(r) emitted by the accretion disk, and hence obtain the effective temperature 
profile T e g(r) = [_F(r)/er] 1//4 , where a is the Stefan-Boltzmann constant. If the disk emits as a 
true blackbody at each radius, it is a simple matter to calculate the total spectral luminosity 
L v dv. By comparing this quantity with the spectral flux F v dv received at Earth, one obtains 
an estimate of rf n cosi/d 2 (essentially the projected solid angle of the disk), where i is the 

The inner edge may not be very pronounced; compare, e.g., the truncated disk model |40j 
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inclination angle and d is the distance to the source. In a few BH binaries, sufficiently reliable 
estimates of i, d and M are available, and thus an estimate of r in is obtained. Identifying 
r- m with r^, one then obtains a. 

A major weakness of this method is that a number of effects will cause the spectrum of an 
accretion disk to deviate from a blackbody. Besides, the method requires accurate estimates 
of M, i, and d. Therefore, spin estimates obtained by this method should be treated with 
caution. 

Quasiperiodic Oscillations 

For some BH binaries, the power spectrum of intensity variations shows one or two 
peaks (more like bumps in some cases) at frequencies of a few hundred Hz. The peaks are 
relatively broad, indicating that they do not correspond to coherent oscillations but rather 
to quasiperiodic oscillations (QPOs). 

One possibility is that the QPO with the highest frequency in each BH binary corresponds 
to the circular Keplerian frequency of gas blobs at some characteristic radius; it is plausible 
that this radius corresponds to the inner edge of the disk. Using equations (I4.16P and (I4.19p . 
one can express the Keplerian frequency of the ISCO as 

^is = -jjF{a/M), (5-14) 

where F(x) is a known function of x — a/M. Assuming that r- m = r; s , one can use this 
method to estimate a provided an estimate of M is available. The method has been applied 
to a few BH binaries (for references, see [8]). Recently, there has been tentative evidence 
for QPOs with a period of 17 minutes in the infrared emission (Genzel et al. |41| ) from Sgr 
A*, the supermassive BH in the Galactic Center . If the QPOs correspond to the Keplerian 
frequency at any radius r > r- m , then the BH must be rotating with a > 0.5 [8J. A number 
of QPO frequencies in the X-ray flares from Sgr A* have been identified by Aschenbach et 
al. [12]. Their analysis reveals that the emission from the inner parts of the accretion disk 
is quite close to the BH horizon and they find a ~ 0.99. 

Relativistic Iron line 

A strong broad spectral line in the X-ray spectrum of the active galactic nucleus (AGN) 
MCG-6-30-15 was recently discovered [43J. They interpreted the line as fluorescent iron Ka 
emission from cool gas in the accretion disk. Similar broad lines were seen in a few other 
AGN and X-ray binaries. Whereas the rest energy of the iron line is 6.4 keV, the observed 
line extends from about 4 to 7 keV. This broadening is due to Doppler blue- and red-shifts 
as well as to the gravitational redshift. 

The line width and its shape, among other factors, depend on the radius range over 
which the emission occurs and in particular on the position of the innermost radius of the 
disk which in turn depends on a (since r in = r is ). 

Given a system with a broad iron line, and assuming that the radiating gas follows 
Keplerian orbits with radii r > r- ls , one can fit the shape of the line profile by adjusting a, 
i, and the emissivity function; the latter is usually modeled as a power law in radius, r _/3 , 
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motivated by the standard disk model [H]. The effect of a is particularly dramatic. As the 
BH spin increases, the inner edge of the disk comes closer to the horizon and the velocity of 
the gas increases substantially. This gives a wider range of Doppler shifts, as well as a larger 
gravitational redshift. The detection of such extreme levels of broadening may be taken as 
a strong indication of a rapidly spinning BH. 

In the case of MCG-6-30-15, the data confirm that the emission comes from a relativistic 
disk and at least some of the data sets can be interpreted in terms of a rapidly spinning 
BH. Assuming that there is no emission from within the ISCO, Reynolds et al. [15] estimate 
a > 0.93. Among BH binaries, the source GX 339-4 shows a broad iron line which seems to 
indicate a > 0.8 (Miller et al. 05]). 

In spite of some weaknesses, the method has the advantage that it requires no knowledge 
of the BH mass or of the distance, and it solves for the disk inclination % using the same line 
data from which a is estimated. 

The variability of the line with time means that it will be challenging to make fundamental 
tests of gravity with this method. On the other hand, the variability could provide interesting 
opportunities to study disk dynamics and turbulence [HI HH] (movies courtesy P. Armitage). 

6 Supermassive Black Holes 

It is now widely accepted that quasars and active galactic nuclei are powered by accretion 
onto massive black holes [^9l [50] [51] . Further, over the last few years there has been increasing 
evidence that massive dark objects may reside at the centers of most, if not all, galaxies 
[521 [53]. I n several cases, the best explanation for the nature of these objects is that they are 
supermassive black holes, with masses ranging from 10 6 to 10 10 solar masses. Comprehensive 
lists of about 30 supermassive BHs at galactic centers may be found in [2J EU 1551 15T)] . 

6.1 Masses and Radii 

The main criterion for finding candidates for such black holes is the presence of a large mass 
within a small region. The mass and the size are estimated using mostly the following three 
methods: gas spectroscopy, maser interferometry, and measuring the motion of stars orbiting 
around the galactic nucleus. 

6.1.1 Gas Spectroscopy 

An example of measurements via gas spectrography is given by the analysis of the Hubble 
Space Telescope observations of the radio galaxy M 87 [57J [5S] . 

A spectral analysis shows the presence of a disklike structure of ionized gas in the inner- 
most few arc seconds in the vicinity of the nucleus of M 87. The velocity of the gas measured 
by spectroscopy at a distance from the nucleus of the order 20 pc ~ 6 x 10 14 km, shows that 
the gas recedes from us on one side, and approaches us on the other, with a velocity differ- 
ence of about 920 km s" 1 . This leads to a mass of the central object of ~ 3 x 10 9 M , and 
no form of matter can occupy such a small region except for a black hole. This is the most 
massive black hole ever observed. 
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6.1.2 Maser Interferometry 



A clear and compelling evidence for black holes has recently been discovered in the radio 
regime: H2O masers orbiting compact supermassive central objects. The structure of ac- 
creting material around the nearby galaxy NGC 4258 (d ~ 6.4 Mpc) has been studied in 
detail [59j [60j [61] . with the aid of very long baseline interferometry (VLBI), which provides 
an angular resolution as fine as 200 /ias (microarcseconds) at a wavelength of 1.3 cm and a 
spectral resolution of 0.1 kms -1 or less, radio interferometry measurements have shown that 
the gas follows circular orbits with a nearly perfect Keplerian velocity profile (v oc r -1 / 2 , 
see Fig. 1 in [8j). Furthermore, the acceleration of the gas has been measured and it too 
is consistent with Keplerian dynamics (Bragg et al. |62j). From these measurements it is 
inferred that there is a dark object with a mass of 3.5 x 1O 7 M confined within ~ 0.13 pc 
= 4 x 10 12 km of the center of NGC 4258. The case for this dark mass being a BH is quite 
strong. 



6.1.3 Virial Mass 

A simple method to obtain a mass estimate is based on the virial theorem. It uses the 
measured velocity dispersion of stars in the central region (ref. [Bj, p. 1007). For simplicity, 
we restrict attention to a spherical cluster of radius R with N stars, each of mass m, so the 
total mass of the bulge is M = Nm. 

The time-averaged kinetic and potential energies of stars in the galaxy's central region 
are related by the equation (ref. [6J, p. 54-56) 

where I is the region's moment of inertia. If the galaxy is in equilibrium, then (d 2 I/dt 2 ) = 0, 
resulting in the usual statement of the virial theorem 

-2(K) = (U) . (6.2) 

Furthermore, for a large number of stars, the central bulge will look the same (in statistical 
sense) at any time, and the time averaging can be dropped. So for iV stars of equal mass, 
we find 

The average velocity squared is 

- ^E^ = V) = «> + «> + <«5> = 3 W = 3 ^ > ( 6 - 4 ) 
where a r is the dispersion in the radial velocity. 

Velocity dispersion 

For a region in space with a large number of stars N, a x measures the spread in the 
x-component of the peculiar velocities of stars and is defined as 

1/2 



1 

x N 



E 



(6.5) 



28 



It is equal to the standard deviation of the velocity distribution in the special case 
when (v x ) = 0. 



Using the (approximate) gravitational potential energy of a spherical distribution of the total 
mass (exercise) 

3 GM 2 



equations (16. 2p and (16.31) yield 



u = ~i—' < 6 - 6 > 



M virial = H£jfr , (6.7) 



where the mass obtained in this way is called the virial mass. 

This equation can be used to estimate a virial mass for the central BH of M31 (An- 
dromeda). The central radial velocity dispersion is measured to be approximately 240 km/s 
within 0.2 as (arcseconds). Given the distance to Andromeda of 770 kpc, 0.2 as corresponds 
to R ~ 0.8 pc. This gives a total mass of roughly 

M virial = 6 x 1O 7 M , (6.8) 

within a sphere of radius 0.8 pc. 

This is, of course, just an order of magnitude estimate. Other estimates for the mass of 
the supermassive BH of M31 range from about 10 6 to 1O 7 M . 



6.2 Sagittarius A* 

Perhaps the most convincing evidence for a black hole comes from the center of our own 
galaxy that coincides with the enigmatic strong radiosource Sgr A*. The existence of a dark 
massive object at the center of the Galaxy has been inferred from the motions of stars and 
gas in its vicinity. The motions of the stars were observed and recorded for many years by 
two independent groups [531 EH ESI ESI EZ] • High- resolution infrared observations made it 
possible to follow the orbits of individual stars around this object (Schodel et al. [63, 64J; 
Ghez et al. [661 EZ])- The movies (1st movie courtesy R. Genzel; 2nd movie courtesy A. 
Ghez [68]) show time-elapsed images of the Galactic Center region revealing the (eccentric) 
orbits of several stars. 

The projected positions of the star S2(S0-2) that was observed during the last decade 
[BUI El] suggest that S2(S0-2) is moving on a Keplerian orbit with a period P of 15.2 yr 
around Sgr A* and the estimated semimajor axis a of the order of 4.62 mpc. 

Then, neglecting the star mass, Kepler's third law (15.21) gives 

47r 2 a^ 

M SgrA * = = 3 ' 7 x 106M © • ( 6 - 9 ) 

The salient feature of the new adaptive optics data is that, between April and May 2002, 
S2(S0-2) apparently sped past the point of closest approach with a velocity v ~ 6000 km/s at 
a distance of about 17 light-hours [63] or 123 AU from Sgr A*. This implies that an enormous 
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mass of the central object is concentrated in a very small volume, strongly suggesting that 
the object is a BH. 

Another star, SO-16 (S14), which was observed during the last few years by Ghez et al. 
[67] with the Keck telescope in Hawaii, recently made a spectacular U-turn, crossing the 
point of closest approach at an even smaller distance of 8.32 light-hours or 60 AU from Sgr 
A* with a velocity v ~ 9000 km/s. Ghez et al. thus conclude that the gravitational potential 
around Sgr A* has an approximately r _1 form, for radii larger than 60 AU, corresponding 
to 1169 .Rsch, where .Rsch = 2M = 0.051 AU for M = 2.6 x 1O 6 M . 

6.3 Supermassive BHs in Active Galactic Nuclei 

Many galaxies possess extremely luminous central regions, with luminosity (in particular 
the luminosity of X-ray radiation) exceeding the luminosity of ordinary galaxies by several 
orders of magnitude. These luminous central regions are called active galactic nuclei (AGN). 
To this class belong the so-called quasistellar objects (QSO) and <yuasars[§. 

What makes these galaxies "active" is the emission of enormous amounts of energy from 
their nuclei. Moreover, the luminosities of active galactic nuclei fluctuate on very short time 
scales - within days or sometimes even minutes. The time variation sets an upper limit 
to the size of the emitting region. For this reason, we know that the emitting regions of 
active galactic nuclei are only light-minutes or light-days across, making them less than one 
ten-millionth the size of the galaxy in which they sit. 

How could a luminosity hundreds of times that of an entire galaxy be emitted from a 
volume billions of times smaller? Of all proposed explanations, only one has survived close 
scrutiny: the release of gravitational energy by matter falling towards a black hole [49, 69J. 
Even using an energy source as efficient as gravity, the black holes in active galactic nuclei 
would need to be supermassive in order to produce the luminosities of quasars. 

6.3.1 Eddington Limit 

The most efficient way of generating energy is by the release of gravitational potential energy 
through mass accretion. For example, a simple calculation shows that for matter falling 
straight down onto the surface of a 1.4 M & neutron star, about 21% of the rest mass is 
released. This is almost 30 times larger than the energy that hydrogen fusion can provide. 

The efficiency may be even larger if the matter is accreted through an accretion disk of 
a rotating black hole. The accretion luminosity (i.e., the radiation energy released per unit 
of time) generated by a mass accretion rate, M, through the disk may be written as 

W = r)M, (6.10) 

where rj is the radiative efficiency of the disc equal to the binding energy of the innermost 
stable circular orbit per unit rest mass. As we have shown in section 14.1.31 0.0572 < 77 < 
0.423, with the lower and upper bounds for a nonrotating and for a maximally rotating BH, 
respectively. 

6 Quasar is short for quasistellar radio source. Quasars are radio-loud whereas QSOs are radio-quiet . 
These names are sometimes confused in the literature. Both terms, QSO and quasar, are often used to refer 
to both types of objects [5J. 
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However, the radiation is interacting with the accreting gas and there is a limit to the lu- 
minosity above which the radiation pressure, acting against gravitational attraction, exceeds 
gravity and thereby stops the accretion. 

Consider a fully ionized hydrogen plasma accreting near the surface of a compact object 
of mass M. The upward force on the infalling matter is mainly due to the interaction of 
radiating photons with electrons in the plasma. If the photon luminosity is L, the number 
of photons iV 7 crossing unit area per unit time at radius r is 

L 

e r 47rr 2 



N i = -TZT^ ( 6 - n ) 



where e r is the mean energy transferred radially per collision. The number of collisions per 
electron per unit time is aN^, where a is the photon-electron scattering cross section. The 
force per electron is just the rate at which the momentum is deposited radially per unit time, 
so we multiply by p r = e r to obtain 

F = aN lPr = ^. (6.12) 
In order for accretion to occur, the gravitational force per electron (acting via the proton) 

*i- = < 613 > 

must exceed the radiation force (16.121) . Here, M. ~ M is the enclosed mass at radius r. 
Equating the radiation force (16.121) with the gravitational force (16.131) sets an upper bound 
to accretion luminosity known as the Eddington limit 

Lem = 4^. (6.14) 
a 

The dominant electron-photon process in a highly ionized hydrogen gas is the scattering of 
photons by free electrons ( Thomson scattering) with the cross section 

a T = — (—} = 6.65 x 1(T 25 cm 2 . (6.15) 



3 V m ^, 
This gives 

L Edd = 1.3 x 10 38 (J^-) erg s" 1 . (6.16) 
Relativistic Eddington Limit 

In a strong gravitational field the expression (16 . 1 3[) for the gravitational force is no longer 
valid. Here we derive the general relativistic expression for the force per electron (acting via 
the proton) under the assumption that the energy density of the fluid is still nonrelativistic, 
i.e., p = m p n and p <^ p. 

Consider a spherical shell of radius r with thickness dr. The force acting on the shell 
element of surface area dS is 

dF = dSdp = dV{l-2M/r) l/2< ^- } (6.17) 
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where 

dV = dr dS(l - 2M /r)- 1/2 (6.18) 

is the proper volume element and dp is the difference between radial pressures at r and r + dr. 
Next, we substitute dp/dr by the right-hand side of the TOV equation (I2.18P in which we 
neglect the pressure (since p C p by assumption). According to ( 12. lip , the number of 
particles in the shell element is dM = ndV. Then, the gravitational force per particle is 

F^ = ^ = -=^l(l- 2 M/r)-V. (6.19) 

Equating |-Fgr av | thus obtained with the photon force ( 16.121) . we find 

Airm v M 



>Edd 



a 



(1 -2M/r)- 1 ' 2 . (6.20) 



Hence, compared with ( 16.141) . the Eddington limit in the strong gravitational field is larger 
by a relativistic correction factor. 

6.3.2 Radius and Mass Estimate 

The luminosity of a typical quasar varies in time with a typical period of 1 hr. This sets the 
upper limit to the radius of about 

R ~ 7AU = 1.1 x 10 9 km. (6.21) 

Considering that AGNs are the most luminous objects known, this is an incredibly small 
size. 

The typical quasar luminosity of 5 xlO 46 erg s _1 is equivalent to more than 500 galaxies 
of the Milky Way size! Now, the constraint that the luminosity must be less than the 
Eddington limit, L < L-^dd, provides a lower limit for the mass of the central object 

M > Mo = 3.3 x 1O 8 M . (6.22) 

The mass limit M$ is quite close to the mass Mbh = R/{2G) = 3.7 x 1O 8 M of a BH the 
Schwarzschild radius of which is equal to the radius of the quasar R estimated above. This 
fact supports the idea that supermassive BHs are responsible for powering AGNs. (Carroll 
and Ostlie [6], p. 1181) 

The estimated value of the lower mass limit is somewhat smaller if the general relativistic 
correction to the Eddington limit is taken into account. Using the relativistic expression 
( I6"20|) . we find 

M>M »l^ 1 + 4i-2^)- 214Xl ° 8M - < 6 - 23 > 



32 



6.3.3 Estimate of the AGN Efficiency 



According to equation ( 16.101) . the radiative efficiency r\ of an accretion disk is defined as the 
energy it radiates per unit accreted mass. As shown in section 14.1.31 77 equals the binding 
energy of gas at the ISCO, which in turn depends on the BH spin parameter a. 

In a typical accretion system, one can easily measure the accretion luminosity L^ sc (pro- 
vided the distance is known), but one practically never has an accurate estimate of the mass 
accretion rate M, so one cannot calculate rj for an individual AGN with the precision needed 
to estimate a. However, a rough estimate can be made for an average AGN. From obser- 
vations of high redshift AGN, one can estimate the mean energy radiated by supermassive 
BHs per unit volume of the universe. Similarly, by taking a sample of supermassive BHs in 
nearby galaxies, one can estimate the mean mass in BHs per unit volume of the universe at 
present. Assuming that supermassive BHs acquire most of their mass via accretion (a not 
unreasonable hypothesis), one can divide the two quantities to obtain the mean radiative 
efficiency of AGN. The current data suggest an efficiency 77 ~ 0.1 — 0.15 for supermassive 
BHs on average (Elvis et al [70], Yu and Tremaine [71]). Such large values of rj are possible 
only if supermassive BHs have significant rotation. 

It should be noted that this is only a statistical result for the population of AGN as a 
whole, and the method does not say anything about the rotation of any specific BH. 

7 Intermediate Mass BHs 

So far we have been concerned with either the stellar mass BHs (M ~ 3.5 — 2OM ) or the 
supermassive BHs (M > 1O 6 M ). Are there BHs in the intermediate mass range 10 2 — 
1O 4 M ? There is a priori no reason for such BHs not to exist. 

There is a tentative evidence based on the Eddington limit that the intermediate mass 
BHs do exist. In several nearby galaxies a number of ultraluminous X-ray sources have been 
detected [72] • As their luminosities of the order of 10 40 or more exceed the Eddington limit 
(I6.16P of 10 M Q BHs, these objects are argued to be good candidates for the intermediate 
mass BHs [75]. However, large apparent luminosity may, under special conditions, be pro- 
duced by stellar mass BHs, e.g, by a supercritical accretion (for details and references, see 
Poutanen et al. [73]). Hence, it is at present unclear what exactly the ultraluminous X-ray 
sources are. Dynamical mass measurements would settle the issue but, unfortunately, none 
of the sources has a binary companion to provide a robust mass estimate. 

8 Observational Evidence for the Horizon 

To prove that a BH-like object is indeed a BH, one needs to demonstrate that it possesses 
an event horizon. The major tests for a BH horizon are based on i) advection-dominated 
accretion flows (ADAF), ii) X-ray bursts, and iii) direct imaging [51 175] 176] 
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8.1 ADAF 



Advection dominated accretion flows [77J, [78] describe accretion with very low radiative 
efficiency in which the energy released by viscosity friction removing the angular momentum 
from the accreting matter is not radiated away but stored in the flow. If an ADAF forms 
around a BH, the stored energy will be lost forever under the event horizon, whereas if the 
accreting object is a NS, this energy must be radiated away once matter lands on its surface. 
Therefore, the argument goes, BHs should be dimmer than NSs if an ADAF is present in 
both cases. 

Narayan, Garcia, and McClintock [79J suggested that precisely such a comparison could 
be done using X-ray binaries. Most of the known stellar-mass BH candidates are in a class 
of X-ray binaries called X-ray novae. These systems are characterized by a variable mass 
accretion rate, and tend to spend most of their time in a quiescent state with a very low 
accretion rate M and accretion luminosity L acc . Only occasionally do they go into outburst, 
when they accrete with high M and become bright. Spectral observations of quiescent BH 
binaries can be explained in terms of an ADAF. Narayan et al. [79] compared quiescent 
luminosities of X-ray binaries supposed to contain BHs with those of neutron-star X-ray 
binaries and realized that, in accordance with the prediction of the ADAF model, systems 
containing black-hole "candidates" are dimmer. They came to the conclusion that they found 
evidence for the presence of event horizons. However, this conclusion has been challenged 
(for details and references, see [51 175]). 

8.2 X-ray Bursts 

In some X-ray binaries, X-ray bursts are observed in addition to the quiescent X-ray lu- 
minosity. In a typical X-ray burst luminosity increases up to nearly the Eddington limit 
in less than a second, and the flux then declines over a period of a few seconds or a few 
tens of seconds. Remarkably, no X-ray bursts of this kind were observed in any BH binary. 
Narayan and Heyl [80] argued that the lack of bursts is a strong evidence for the horizon in 
BH candidates. 

The explanation for such an absence of X-ray bursts is quite simple. At present, it is 
widely accepted that X-ray bursts arise from thermonuclear detonation of accreted material 
[801 18T] . Accreted material builds up on the surface of the compact object and is com- 
pressed by the object's gravity. After sufficient material accumulates, it undergoes unstable 
thermonuclear burning, which we observe as an X-ray burst. 

A key point is that the compact object must have a surface. Material cannot accumulate 
on an event horizon, and so no bursts can come from an X-ray binary whose compact object 
is a BH. For more details, references, and critical discussion, see [H [75] 

8.3 Direct Imaging 

The most promising line of search for a direct evidence is to construct an image of the region 
near the event horizon using interferometry. 

For defmiteness, consider a nonrotating BH with a horizon at radius R. Because of 
strong gravitational lensing in the vicinity of the BH, a distant observer will see an apparent 
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boundary of the BH at a radius of ?>\f?>R/2 (82]. Rays with impact parameters inside this 
boundary intersect the horizon, while rays outside the boundary miss the horizon. The 
angular size of the boundary, e.g., for Sgr A* with M ~ 4 x 1O 6 M and at a distance of 8 
kpc, is ~ 0.02 mas which is not beyond reach. The supermassive BH in the nucleus of the 
nearby (15 Mpc) giant elliptical galaxy M87, with a mass of 3 x 1O 9 M and an expected 
angular size of ~ 0.01 mas, is another object of interest |83j. 

The best angular resolution achievable today is with radio interferometry, where angles 
less than 1 mas are routinely resolved. In the not too distant future, it should be possible to 
operate interferometers at wavelengths A < 1 mm and with baselines as large as the diameter 
of the Earth b ~ 10 4 km. For details, see Falcke et al. [82] . 

9 Alternatives to Supermassive BHs 

Given the accumulated evidence for supermassive compact objects ranging from a few 10 6 M e 
to a few 1O 9 M , the existence of black-hole-like objects is beyond doubt J53j[56]. What still 
remains an issue is whether these supermassive objects are BHs with the Schwarzschild (or 
Kerr) metric describing the physics of the interior, or some other objects built out of more 
or less exotic substance but with a regular behavior in the interior. A standard astrophysical 
scenario in the form of a compact cluster of dark stars (e.g., neutron stars or brown dwarfs) 
although not entirely excluded, is quite unlikely. It has been demonstrated that, in the 
case of NGC 4258 and our Galaxy, such a cluster would be short-lived and would either 
"evaporate" or become a BH in much less time than the age of the Galaxy [84J. 

A number of alternatives to classical BHs have been proposed with no singularities in the 
interior. The three representative models are: 1) Neutrino (or neutralino) stars, 2) Boson 
stars, 3) Dark energy stars 

9.1 Neutrino Stars 

Here we use the term neutrino stars as a generic name for any degenerate fermion star com- 
posed of neutral weakly interacting fermions, e.g., neutrinos or supersymmetric fermionic 
partners such as neutralinos, gravitinos, and axinos. The simplest model proposed for su- 
permassive compact objects at the galactic centers is a self-gravitating degenerate fermion 
gas composed of, e.g., heavy sterile neutrinos [85l IBS EH ES ES]. Sterile neutrinos in the 
keV mass range have recently been extensively discussed as dark-matter candidates [901 ED] 

As we have seen in the example of a neutron star, a self-gravitating ball of degenerate 
fermionic matter is supported against gravitational collapse by the degeneracy pressure of 
fermions due to the Pauli exclusion principle, provided the total mass is below the OV limit 
(I3.14p . with an m~ 2 functional dependence on the fermion mass m. 

Let us assume that the most massive objects are sterile neutrino stars at the OV limit 
with Mqv — 3 x 10 9 Mq, such as the supermassive compact dark object at the center of the 
radio-galaxy M87 [57]. Using (I3.14p we find that the neutrino mass required for this scenario 
is 

m ~ 14 keV for g = 2, 

m ~ 12keV for g = 4. (9.1) 
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From (13.141) and (13. 151) it follows that a neutrino star of mass Mov = 3 x 10 9 M & would have 
a radius R y = 4.4466i? S ch = 3.9396 x 10 10 km = 1.52 light-days, where R Sch = 2GM ov is 
the Schwarzschild radius for Mov- Thus, at a distance of a few Schwarzschild radii away 
from the supermassive object, there is little difference between a neutrino star at the OV 
limit and a BH, in particular since the last stable orbit around a BH already has a radius of 

3 -Rsch- 

Now, we wish to investigate the possibility that this scenario can be extrapolated to 
explain all observed supermassive BH, and in particular Sgr A*, the supermassive compact 
dark object at the Galactic center. As the mass of Sgr A* is by about a factor of 10 3 lower 
than the OV limit above, we can use the Newtonian mass-radius scaling relation ( 12.321) 
derived from the Lane-Emden equation f!2.28j) with the polytropic index n = 3/2. We find 

f!y /3 . 7,40 6 xl0^^X) 8/ 7i) 2/ 7^) ,/3 AU. (9.2) 
m s/3 M i/3 y g J \ m J \gj \ M J y ' 

Here, 1 mpc = 206.265 AU. The degeneracy factor g = 2 describes either spin 1/2 fermions 
(without antifermions) or spin 1/2 Majorana fermions. For Dirac fermions (or spin 3/2 
Majorana fermions), we have g = 4. Using the canonical value M = 2.6 x 1O 6 M , we find 
R ~ 5200 AU to be compared with the observational upper limit |67J R < 60 AU for Sgr A*. 
Hence, a neutrino star made of m = 12 — 15 keV neutrinos is ruled out as an explanation 
for Sgr A*. 

To be able to explain Sgr A* as a neutrino star, we need a minimal fermion mass m m i n 
= 63.9 keV/c 2 for g = 4. The maximal mass of a neutrino star made of these fermions, as 
given by the OV limit (l3TT4"j) . is [92] 

M™ x = 1.083 x 1O 8 M . (9.3) 

Clearly, a neutrino star scenario which would cover the whole mass-radius range of com- 
pact supermassive galactic centers is ruled out. However, a "hybrid" scenario [H2] in which 
all supermassive compact dark objects with masses M > Moy x are black holes, while those 
with M < Mqy x are neutrino stars, is not excluded. 

At first sight, such a hybrid scenario does not seem to be particularly appealing. However, 
it is important to note that a similar scenario is actually realized in nature, with the co- 
existence of neutron stars with masses M ^$ 3M Q , and stellar- mass BHs with mass M > 3M , 
as observed in binary systems in the Galaxy. 

An indirect support for this scenario is the absence of clear evidence for intermediate 
mass BH candidates, which is difficult to explain in the conventional BH scenario (in which 
BHs are all baryonic). If the hybrid scenario were realized in nature, the intermediate mass 
neutrino stars would exist but having very large radii would be rather dilute and hence 
difficult to detect. 

9.2 Boson Stars 

Boson stars are static configurations of self-gravitating (with or without self-interaction) 
complex scalar fields. In self-interacting scalar field theories, such as 4 theory, there are 
cases where a homogeneous condensate is a stable ground state, known as the Bose-Einstein 
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(BE) condensate. Hence, a boson star being basically a self-gravitating BE condensate may 
be regarded as a T — > limit of a self-gravitating boson gas [93]. The ground state of a 
condensed cloud of charged bosons of mass m, interacting only gravitationally and having 
a total mass M below a certain limit of the order Mjj/m, is a stable spherically symmetric 
configuration [M] which is usually referred to as mini-soliton star [95j or mini-boson star 
[96] . For a recent review and a comprehensive list of references, see Schunck and Mielke [97J. 

The gravitational collapse of a self-interacting BE condensate is prevented by a repulsive 
self-interaction, e.g., in the form of A|$| 4 . That makes it astrophysically interesting as its 
maximal mass is ~ mj^/mfj, hence comparable with the mass of a neutron star or a neutrino 
star. 

Seidel and Suen have demonstrated that boson stars may be formed through a dissipa- 
tionless mechanism, called gravitational cooling [98J. Boson stars have recently attracted 
some attention as they may well be candidates for nonbaryonic dark matter [99] . 

The action of the gravitationally coupled complex scalar field $ reads 

Sbs = j d 4 x v^det^ (-^- + £ KG ) , (9.4) 

with the Klein-Gordon Lagrangian 

Ckg = <r(fy$*)($,$) - m 2 |$| 2 - Um 2 ) . (9.5) 
By varying the action with respect to $ one obtains the Klein-Gordon equation 

and the variation with respect to g^ v yields Einstein's equations flA.23j) in which 



T„ v = {d^*){d v <5>) + {d^){d u $*) - g^C KG , . (9.7) 

Although this energy-momentum tensor is not in the perfect fluid form fIB.ll we can still 
identify the radial pressure p and the density p with — T r r and Tq , respectively. 

Given the interaction potential U, one has to solve a coupled system of Klein-Gordon 
and Einstein equations in spherically symmetric static spacetime. Unlike in the case of a 
fermion star, where the information about matter was provided by the equation of state, the 
matter here is described by the KG equation (I9.6p . The stationarity ansatz 

$(r,t) = -^(r) e -^ (9.8) 

describes a spherically symmetric bound state of the scalar field, where the frequency u is 
determined by the asymptotic condition ip(r) — >• as r — > oo. 

It turns out that even the simplest case U — has a nontrivial ground-state solution, 
called a mini-boson star. The gravitational collapse of a mini-boson star is prevented by 
Heisenberg's uncertainty principle. This provides us also with crude mass estimates: For 
a boson to be confined within the star of radius Ro, the Compton wavelength has to sat- 
isfy A$ = 2-K/m ~ 2Rq. On the other hand, the star's radius should be larger than the 
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Schwarzschild radius, R$ > i?sch = 2M/mpi in order to avoid instability against complete 
gravitational collapse. In this way we obtain an estimate 

MK^mlJm. (9.9) 

This upper bound is slightly larger than the so-called Kaup limit, 

M max = 0.633m P1 /m = 8.4639 x KT 10 \^~\ M Q ( 9 - 10 ) 

obtained numerically [96]. Clearly, mini-boson stars are irrelevant to astrophysical context, 
unless the boson mass m is ridiculously small. 

This result was later extended by Colpi et al. |10U] for interacting bosons with a quartic 
self-interaction 

[7 = A|$| 4 . (9.11) 
In this case, the maximal mass is obtained as [10(3, HOlj 

M -^^^ l H^f) M ° (9 - 12) 

quite similar to the fermion star mass limit. The effective radius of an interacting boson 
star in the Newtonian regime turns out to be independent of central density and in a good 
approximation is equal to 

R cS = -^V\^ = 1.5096VA ( 2 km. (9.13) 

Hence, boson stars with a quartic self-interaction all have the same radius. This behavior 
is quite different from the fermion star case, where the radius scales with the star mass as 
R ~ M 1 / 3 (see equation (IO> ). 

The results are summarized in the table 



Object 


Maximum 

Mass 


Effective 
Radius 


Fermion Star 
Mini-Boson Star 
Boson Star 


M ov = 0.384mp L /m 2 
M Kaup = 0.633m PL /m 
M max = l/(v / 87r3) v ^m| 1 /m 2 


Row — 3.36m P1 /m 2 
1/m 

R e s = 7r/(-\/87r)\/Ampi/m 2 



Could boson stars fit the whole range of compact supermassive galactic centers? To 
answer this question, we proceed as in the case of neutrino stars. We assume that the 
most massive objects are boson stars with maximal mass, i.e., M max ~3x 10 9 M Q , for the 
supermassive compact dark object at the center M87. From this we find 




(9.14) 
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which also fixes the radius given by (19.131) 



R = 4.53 x 10 10 km = 303 AU , (9.15) 

quite close to the radius Roy of the neutrino star discussed in section 19.11 Since the radius 
of a boson star does not depend on its mass, this value should also fit the radii of all 
supermassive BHs. However, it obviously does not fit the observed radius limit of Sgr A*, 
R < 60 AU, and even less it fits the radius of a typical AGN of 7 AU (section r6.3.2j) . 

Hence, boson stars with a quartic self-interaction are ruled out as a unique explanation 
for compact supermassive objects at galactic centers and AGN. Again, it is not excluded 
that some of the galactic centers harbor a boson star. 

9.3 Dark-Energy Stars 

Dark energy (DE) is a substance of negative pressure needed for an accelerated cosmology. 
In order to achieve acceleration, DE must satisfy an equation of state p^ c = wp^e, where w 
may depend on the cosmological scale a but at present (a = 1) w < — 1/3. In addition, one 
also expects that DE satisfies the dominant energy condition, in which case the equation 
of state satisfies w > — 1. Nevertheless, matter with w < — 1 dubbed "phantom energy" 
|102[ I103j has recently received attention as it might fit the most recent SN la data slightly 
better than the usual DE. 

In cosmological context, pde is homogeneous and it is normally assumed that DE does 
not cluster. However, it is not excluded that owing to gravitational instability, small inho- 
mogeneities (analogous to dark- matter inhomogeneities) grow and build structure, e.g., in 
the form of spherical configurations. Hence, we define dark-energy stars as spherical solu- 
tions to Einstein's equation with matter described by the DE equation of state [104[ [105J. 
It is worth noting that DE could be described by a scalar field theory (quintessence) with a 
suitably chosen interaction potential and/or with a noncanonical kinetic energy term in the 
Lagrangian. This scalar field theory would correspond to an effective equation of state with 
desirable properties |106j . In this way, DE stars could be regarded as boson stars but with 
a rather unusual self-interaction. 

Here, we briefly discuss two examples of DE stars: de Sitter gravastars and Chaplygin 
gravastars. 

9.3.1 De Sitter Gravastars 

The simplest example of DE star is a gravitating vacuum star or a gravastar. Chapline et al. 
|1071 1108] put forth an interesting proposal based on analogies to condensed matter systems 
where the effective general relativity was an emergent phenomenon. Specifically, assuming 

1 /2 

the Schwarzschild exterior, they suggested that the sphere where the lapse function £ = g ' 
vanished marked a quantum phase transition, £ increasing again at r < 2GM. As this 
required negative pressure, the authors of [108] assumed the interior vacuum condensate to 
be described by de Sitter space with the equation of state p = —p. 

Subsequently, the idea of gravitational vacuum condensate, or 'gravastar', was taken up 
by Mazur and Mottola [109[ 1110] and Dymnikova |lllj , replacing the horizon with a shell of 
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stiff matter astride the surface at r = 2GM. Visser and Wiltshire |112j and recently Carter 
[113] also examined the stability of the gravastar using the Israel thin-shell formalism [114]. 

Given the mass M, the interior is described by a solution, similar to the the Einstein-de 
Sitter universe, with constant density p = p up to the de Sitter radius Rds = 2M (actually 
slightly further out). Hence, the interior is de Sitter, p being the vacuum energy density. 
The lapse function is given by 




with 

R*s = \I^Po 1/2 - (9-17) 

In order to join the interior solutions to a Schwarzschild exterior at a spherical boundary 
of radius R, it is necessary to put a thin spherical shell at the boundary with a surface 
density and a surface tension satisfying Israel's junction conditions |114j . As the pressure 
does not vanish at the boundary, it must be compensated by a negative surface tension of 
the shell. 

The gravastar has no horizon or singularity! Its surface is located at a radius just slightly 
larger than the Schwarzschild radius R = 2M + e, e ~ 2lpi- Hence, there is practically 
no observational way to distinguish a Schwarzschild BH from a gravastar. Basically, any 
BH-like object observed in nature may be a gravastar. 

In spite of these attractive features the gravastar has been met with a cool reception. 
Certainly, the assumption of a de Sitter interior presents a quandary: on the one hand, 
the quantum phase transition would suggest that the associated cosmological constant is a 
fundamental parameter; on the other hand, to accommodate the mass range of supermassive 
black-hole candidates, it must vary over some six orders of magnitude. In addition, the 
notorious cosmological constant problem is reversed for gravastars: why is the vacuum energy 
in the interior of a gravastar so much larger than the observed vacuum energy density in 
the universe? If we identify the most massive black-hole candidate observed at the center 
of M87, with mass M max = 3 x 1O 9 M , with the de Sitter gravastar, then p = (9.7keV) 4 , 
to be contrasted with the (10 _3 eV) 4 values wanted for cosmology. Another question is how 
does a gravastar form? The entropy of a stellar mass gravastar is much less than the entropy 
of an ordinary star and this would require an extremely efficient cooling mechanism before 
gravastars could form during stellar collapse [75j . 



9.3.2 Chaplygin Gravastars 

Consider a particular form of DE with a rather peculiar equation of state 

A 

p= . (9.18) 

P 

Equation (19.18p describes the Chaplygin gas which, for p > y/A, has attracted some attention 
as a dark-energy candidate [106[ 1115] . Consider a self-gravitating Chaplygin gas and look 
for static solutions. In particular, we look for static Chaplygin gas configurations in the 
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phantom (w < —1) regime, i.e., when 



p<VA. (9.19) 

We show that these configurations could provide an alternative scenario for compact massive 
objects at galactic centers |116j . Moreover, equation (I9.18P yields the de Sitter gravastar 
solution in the limit when the central density of the static solution approaches the value 
'A. 

Combining (BTIBji with the TOV equations fTZTSl) . and (12191) . one has 

/ _ £\ , ( , 20) 



Anr 2 p. (9.21) 



A J V r(r - 2M) 

dM 

dr 

In Fig. 1 we exhibit the resulting p(r) for selected values of po/^/A. The solutions 
depend essentially on the magnitude of po relative to y/A. In the following we summarize 
the properties of three classes of solutions corresponding to whether p is larger, smaller, or 
equal to \[A. 

i) For po > VA, the density p increases and the lapse function £ decreases with r starting 
from the origin up to the black-hole horizon radius R^, where 2GAi(Rbh) = Rbh- In the 
limit po — > oo, a limiting solution exists with a singular behavior 

7 A 

near the origin. 

ii) For po < v^4, both p and £ decrease with r up to the radius Rq where they vanish. At 
that point the pressure p blows up to — oo owing to ( 19. 18ft . The enclosed mass M. is always 
less than r/(2), never reaching the black-hole horizon, i.e., the radius where 2Ai(r) = r. 

iii) For p = VA, the density p remains constant equal to \/A~ up to the de Sitter radius 
i?ds — 2M. Hence, we recover the de Sitter solution precisely as in the gravastar case. 

As po — > VA from above or from below, solutions i) or ii), respectively, converge to iii) 
except at the endpoint. The lapse function in iii) joins the Schwarzschild solution outside 
continuously, whereas in i) and ii) it happens discontinuously. 

As in the case of a de Sitter gravastar, in order to join our interior solutions to a 
Schwarzschild exterior at a spherical boundary of radius R, it is necessary to put a thin 
spherical shelf] at the boundary with a surface density and a surface tension satisfying Is- 
rael's junction conditions [114J. 

Case ii) , together with iii) , is of particular interest as we would like to interpret the super- 
massive compact dark objects at the galactic centers in terms of phantom energy rather than 
in terms of a classical black hole. It is natural to assume that the most massive such object 
is described by the de Sitter gravastar, i.e., solution iii) (depicted by the solid line in Fig. 0j). 
The radius of this object is Rds equal to the Schwarzschild radius 2M max . Clearly, solutions 



7 It has recently been demonstrated that the joining can be made continuous without the presence of a 
thin shell for a gravastar made of the fluid with an anisotropic pressure |117[ 1118] . 
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Figure 4: Density profile of the Chaplygin gravastar for p /\^A = 1.2 (short dashed), 0.98 
(dotted), 0.9 (long dashed). The limiting singular solution with the r~ 2 / 3 behavior at small 
r is represented by the dot-dashed and the de Sitter gravastar by the solid line. 



belonging to class ii), can fit all masses M < M max . However, for the phenomenology of 
supermassive galactic centers it is important to find, at least approximately, the mass-radius 
relationship for these solutions. This may be done in the low central density approximation, 
i.e., po ^ v^4, which is similar to the Newtonian approximation but, in contrast to the New- 
tonian approximation, one cannot neglect the pressure term in (12.181) . Moreover, as may be 
easily shown, in this approximation M. <C r 3 p, so that the pressure term becomes dominant. 
Next, neglecting 2Ai with respect to r, as in the usual Newtonian approximation, equation 
(I9.20p simplifies to p' = AirAr, with the solution 

P^ofl-^); (9.23) 



which gives a mass-radius relation 

M 16tt 2 



R% 15 



A = constant. (9.24) 



The mass-radius relationship M oc Rq which phantom gravastars obey, offers the prospect 
of unifying the description of all supermassive compact dark objects at the galactic centers, 
as Chaplygin-gas phantom gravastars with masses ranging from M min = 1O 6 M to M max = 
3 x 10 9 M Q . Indeed, assuming that the most massive compact dark object, observed at the 
center of M87, is a Chaplygin phantom gravastar near the de Sitter gravastar limit, with 
-Rmax = 2M max = 8.86 x 10 9 km = 8.21 lhr, the compact dark object at the center of our 
Galaxy, with mass Mqq = 3 x 10 6 M q , would have a radius Rqc = 2.06 lhr if the scaling law 
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(I9.24p holds. This radius is well below the distances of the closest approach to Sgr A* which 
the stars SO-2 (R min = 17 lhr = 123 AU, [631 El]) and SO-16 (R min = 8.32 lhr = 60 AU 
[66| [67] ) recently had and beyond which the Keplerian nature of the gravitational potential 
of Sgr A* is well established. 

A Basic General Relativity 

A.l Geometry 

The geometry of spacetime is described by the metric 

ds 2 = g^dx^dx" , (A.l) 

where g^ v is the metric tensor. The spacetime curvature is defined through the Riemann 
curvature tensor 

RrtT/37 — QraR^p*/ j (A -2) 

defined by 
where 

= \g Xa [d,guX + d u g, x - d x g^\ (A.4) 
are the Christoffel symbols. A contraction of two indices gives the Ricci tensor 

Raf3 = g XV R\av(3 (A. 5) 

and a further contraction of the Ricci tensor gives the scalar curvature 

n = g afi R afi . (A.6) 

A. 2 Covariant Derivative 

A covariant derivative V M of a field Y (in general, Y may be a tensor of any rank) is defined 
by 

BY = (V^dx* = Y.^dx" , (A.7) 

where DY is an infinitesimal difference between the value of the field Y at the point x^ + dx^ 
and the quantity Y{x^) parallelly displaced from x^ to x^ + dx^. Hence, DY consists of two 
parts: one is the change of Y due to the parallel displacement and the other is the difference 
dY = Y(x^ + dx^-) — Y(x^) due to the functional dependence on x M . The latter part is 
basically related to the ordinary partial differentiation. The difference due to the parallel 
displacement is related to the curvature of spacetime and depends on the tensor nature of 
Y. For a scalar, vector, and second-rank tensor one finds: [119] 

• scalar (p 

<P;v = <P,n , (A.8) 
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• vector 

v„„, = - v%y 9 , = v% + v%y . (A.9) 

• tensor 

A -A — V a A — V a A A^ u — A^ u _i_ A au 4- V A^ u ( A 1 0^1 

The covariant d'Alembertian is given by 

□ = <rv„v„ = (v^dit^^) . (A.n) 

Here we have used the usual convention in which a subscript , \x denotes an ordinary 
partial derivative and ; \i denotes the covariant derivative. 



A. 3 Geodesies 

Geodesies are the "shortest possible lines" one can draw in curved geometry. Given a co- 
variant derivative operator V M , we define a geodesic to be a curve whose tangent vector is 
parallel propagated along itself, i.e., a curve whose tangent satisfies the equation 

uV, = 0. (A.12) 

A massive particle moves along a timelike geodesic, w M being its velocity. A massless particle 
moves along a null geodesic in which case the vector w M is null. 



A. 4 Isometries and Killing Vectors 

A vector field k^ 1 that generates one parameter group of isometries, i.e., one parameter group 
of transformations that leave the metric invariant, is called a Killing vector. A Killing vector 
k^ satisfies the Killing equation 

k(p;i/) = kp,;v + ky-^ = 0. (A. 13) 

It is convenient to represent the Killing vector /c M as a differential operator k, 

k = k^. (A. 14) 

For a vector field k, local coordinates can be found such that 

*=A, (A.15) 

where x is one of these coordinates, e.g., x A = x. In such a coordinate system, k^ = 5% and 

*W) = %^ = 0. (A. 16) 

Hence, one can say that k^ is Killing if g^ v is independent of x. 

Examples 

Consider the spherical coordinate system (t,r,9,(f)). Suppose the metric components are 
independent of t and <fi. Then there exist two Killing vectors: 
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• the generator of time translations £ = d/dt 

= ^ = g, , (A.17) 

• the generator of axial isometries ip = 8/ dip. In spherical coordinates 

r = #i\ vv = ^- (A.i8) 

A. 5 Constants of Motion 

The following proposition relates Killing vectors to constants of motion. 

Proposition 1 Let x M be a Killing vector field and let 7 be a geodesic with tangent . Then 
the quantity Xn ufl ^ s constant along 7. 

Proof: We have 

^(XXV = u^Xw + = 0. (A.19) 

The first term vanishes by the Killing equation (1A.13I) and the second term vanishes by the 
geodesic equation (1A.12I) . Hence 

= 0, (A.20) 

as was to be proved. D 
Examples 

Consider a stationary axially symmetric space, as in the example discussed in appendix 
IA.4I with two Killing vectors, £ = d/dt and ip = d/d<ft, corresponding to time translation 
and axial isometries, respectively. In addition, assume the asymptotic flatness, i.e., g^ u — ► 
diag(l, — 1, — r 2 , — r 2 sin 2 8) as r — ► 00. For a particle of mass m moving along a timelike 
geodesic, the two conserved quantities along the geodesies are 

• the energy 

E = mu»Z lt = Z lt <f = g 0li <f=q°\ oo , (A.21) 

• the z-component of the angular momentum 

L = -mu^ = = -g^q" = r 2 sin 2 0g*L . (A.22) 

A. 6 Einstein's Equations 

General relativity relates the geometry of spacetime to matter through Einstein's field equa- 
tions 

~ \g^n = -8tt7V , (A.23) 
where T^ u is the energy-momentum tensor. 
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B Basic Fluid Dynamics 

Consider a perfect gravitating relativistic fluid. We denote by , p, p, n, and a the 
velocity, pressure, energy density, particle number density, and entropy density of the fluid. 
The energy-momentum tensor of a perfect fluid is given by 

Tfiu = (p + p)u^u u - pg^u, (B.l) 

where is the metric tensor with the Lorentzian signature (H ). Hence, in this 

convention, we have 

w% = g^u" = 1. (B.2) 
The particle number conservation is described by the continuity equation 

(Tm>% = -^Ld M (V=# mi") = 0. (B.3) 

The energy-momentum conservation 

= (B.4) 
applied to (IB. II) yields the relativistic generalization of Euler's equation [120j 

(p + p)u u u fl - u -d f ,p + u f ,u u d u p = 0. (B.5) 

B.l Fluid Velocity 

It is convenient to parameterize the four-velocity of the fluid in terms of three-velocity 
components. To do this, we use the projection operator g^ u — t^t u , which projects a vector 
into the subspace orthogonal to the time-translation Killing vector = (1; 0) , where is 
the unit vector ^ 

= — = ^° ■ t = ^ = ^ At0 (B 6) 

We split up the vector in two parts: one parallel with and the other orthogonal to t^. 

= jtfj, + (g^ - t^t u )u v , (B.7) 

where 

7 = • (B.8) 
From (1B.7|) with (1B.6|) we find (see, e.g., appendix of [121] ) 



= 7 



9ojV 3 



«m = 7 ( v^oo! "4= _ u * ) ' ( B - 9 ) 



fl'oo #00 

fa 
. #00 

where 

v i = -y ij v i , v 1 = v i v i , 1 2 = {\-v 2 Y\ (B.10) 
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with the induced three-dimensional spatial metric: 

liS = ^-9ii\ iJ= 1,2,3. (B.ll) 

Since w M and are timelike unit vectors, a consequence of ( IB. 71) is that 7 > 1 and hence 
< v 2 < 1. 

B.2 Hydrostatic Equilibrium 

From Euler's equation (1B.5P we can derive the condition of hydrostatic equilibrium. We can 
use the comoving frame of reference in which the fluid velocity takes the form 

u"=-^-; u„ = -^_. (B.12) 
v^oo v^oo 



In equilibrium the metric is static; all components are independent of time, and the mixed 
components got are zero. Equation (IB. 51) then gives 



(p + p)r>%, = (p + p)^°%Soo = -d^p (B.13) 



1 
2 

or 

-i/2~ 1/2 



fyp = ~(p + p)#oo ^0oo • ( B - 14 ) 



C Basic Thermodynamics 

Consider a nonrotating fluid consisting of iV particles in equilibrium at nonzero temperature. 
A canonical ensemble is subject to the constraint that the number of particles 

J nu>*dE fl = N (C.l) 

should be fixed. The spacelike hypersurface £ that contains the fluid is orthogonal to the 
time-translation Killing vector field In equilibrium is related to the velocity of the 
fluid. 

e = ^; e = (^U 1/2 - (C2) 

It may be shown that those and only those configurations will be in equilibrium for which 
the free energy assumes a minimum |122j . The canonical free energy is defined as [1221 [123]. 

F — M — f 2V£"dE M , (C.3) 

where M is the total mass as measured from infinity. The entropy density is obtained using 
the standard thermodynamic relation 

a= ^(p + p-pn). (C.4) 



47 



C.l Tolman Equations 



The temperature T and the chemical potential [i are metric-dependent local quantities. Their 
spacetime dependence may be derived from the equation of hydrostatic equilibrium ( IB. 141) 
[TOf 1120] and the thermodynamic identity (Gibbs-Duhem relation) 



4 

T 



nd^- — pd—. 
T T 



The crucial condition is that the heat flow and diffusion should vanish [124] 

const, 



E 

T 



(C.5) 



(C.6) 



which may be derived from the physical requirement that the rate of entropy change with 
particle number at fixed energy density should be constant, i.e., 



da 
On 



const, 



(C.7) 



where 'const' is independent of p. From (IC.4I) and flC.5j) we obtain 

da = —dp — T=dn 
T T 



and hence 



Oa 

dn 



E 

T' 



(C.9) 



from which (10.6|) immediately follows. Next, equation (IC.6I) . together with (IB. 141) and (|C5p . 
implies the well-known Tolman equations 



rp 1/2 

T 9oo 



T ; 



1/2 

Moo 



1^0 



(CIO) 



where T and /io are constants equal, respectively, to the temperature and the chemical 
potential at infinity. In a grand-canonical ensemble, T and /io may be chosen arbitrarily. In 
a canonical ensemble, fiQ is an implicit functional of £ because of the constraint (IC.ll) that 
the total number of particles should be fixed. 



C.2 Fermi Distribution 

Consider a gas consisting of iV fermions with the mass m contained within a hypersurface 
S. The equation of state may be represented in a parametric form using the well-known 
momentum integrals over the Fermi distribution function [125] 

r° d 3 q 1 
n ~ 9 J J2vf 1 + e^-/VT ' IS" 11 ) 

f°° d 3 q E 
P ~ 9 J J2vf 1 + e^-/VT ' ( C ' 12 ) 
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p=gT [w? ln(1 + e ~ E/T+ ^ ' ( c - 13 ) 

where E = m 2 + q 2 and T and /x are the local temperature and the local chemical potential, 
respectively, defined by Tolman's equations ( 1C.10I) . By partial integration, the last equation 
may be written as 

f°° d 3 q q 2 1 
P = 9 J (2^3El + e^-^' (C ' 14) 

The integer g denotes the spin degeneracy factor. In most applications we take g = 2 (spin 
up and spin down). Strictly speaking, in each equation (1C.11I) - (1C.14P one should also add the 
antiparticle term which is of the same form as the corresponding right-hand side of (10.110 - 
(1C.14I) with n replaced by — fi. However, the contribution of antiparticles in astrophysical 
objects is almost always negligibl^. 

Equations (1C.4I) and (1C.8I) may be combined to yield another useful thermodynamic 
identity 

dw = Td(-) + -dp, (C.15) 

n n 

with w = (p + p)/n being the specific enthalpy. 



C.3 Isentropic Fluid 

Euler's equation is simplified if one restricts consideration to an isentropic flow. A flow is 
said to be isentropic when the specific entropy a/n is constant, i.e., when 

^(-) = 0. (C.16) 

n 

A flow may in general have a nonvanishing vorticity uj^ defined as 

uj,,u = h^hy [p . a] , (C.17) 

where 

K = 5»-u»u v (C.18) 

is the projection operator which projects an arbitrary vector in spacetime into its component 
in the subspace orthogonal to u^. A flow with vanishing vorticity, i.e., when 

OV = °> ( c - 19 ) 

is said to be irrotational. In the following we assume that the flow is isentropic and irrota- 
tional. 

As a consequence of equation ( 1C.16I) and the thermodynamic identity ( 1C.15I) . equation 
(IB.5I) simplifies to 

u"(wuJ. )V - d„w = 0. (C.20) 

8 One exception is neutrino stars made of Dirac type neutrinos. There, the numbers of neutrinos and 
antineutrinos are equal and separately conserved, hence g=4 
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Furthermore, for an isentropic irrotational flow, equation flC.19j) implies [126] 

(wu^-u - (wuu).^ = 0. (C.21) 

In this case, we may introduce a scalar function ip such that 

wUf, = -d^p, (C.22) 

where the minus sign is chosen for convenience. Obviously, the quantity wu^ in the form 
flC.22j) satisfies equation flC.21j) . Solutions of this form are the relativistic analogue of po- 



tential flow in nonrelativistic fluid dynamics [120J. 
C.4 Degenerate Fermi Gas 

In the limit T — > 0, the Fermi distribution in (1C.11|) - (1C.13I) becomes a step function that 
yields an elementary integral with the upper limit equal to the Fermi momentum 



vV - m 2 = mX. (C.23) 



The equation of state can be expressed in terms of elementary functions of X. With g = 2, 
we find 

riF ^3 i r , 

p = 2 4^ = — m 4 X(2X 2 + l)v / X^TT-ArshX , (C.25) 

Jo { 2rr r otc z l j 



QF q 2 d 3 q 1 



P — 2 / -, 77 = 7:7X1 

There are two important limits: 
• Nonrelativistic limit, I C 1 



X{-X 2 - 1)VX 2 + 1 + ArshX 
3 



(C.26) 



>=^ 4 r + ^- !**■■■)■ (c27) 

Obviously, in this limit, p <^ p. 
Ultrarelativistic limit 

p = ^m 4 (x 4 + X 2 - 1 ln2X...) , (C.29) 

P = 12^ ~ X ' + f ln2X ") • (C ' 30) 

Retaining the dominant terms, these equations yield the well-known equation of state 
for an ultrarelativistic gas 

P=f, (C31) 
which also holds for massless Bose and Fermi gases at finite temperature. 
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C.5 Polytropic Gases 

A gas of particles of mass m is called polytropic if the equation of state may be written in 
the form 

p = Kn T , (C.32) 

where K, is a constant that depends on m. Using this equation and the thermodynamic 
identity flClo]) the energy density p and the entropy density a may also be expressed in 
terms of n. For an isentropic flow, it follows 

p = mn + p~~Y nF - (C.33) 

A degenerate Fermi gas approaches a polytropic equation of state in both the nonrelativistic 
and the extreme-relativistic limits: 

• Nonrelativistic limit. Retaining only the dominant term in ( 1C.28I) we find 

o2/3 4/3 

p = 6 -^^n"\ (C.34) 
5m 

• Ultrarelativistic limit. In this case, the dominant term in (1C.30I) yields 

Ql/3_2/3 

p = 6 -^— n '/\ (C.35) 
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